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ABSTRACT 


The  existing  models  for  treating  radiochemical  fractionation  In 
nuclear  weapon  debris  are  discussed  and  compared.  A  method  which 
extends  the  existing  theory  for  the  case  of  surface  bursts  over 
silicate  soils  Is  developed  and  validated  with  weapons  test  data. 

There  Is  evidence  that  fission  product  absorption  by  soil  and 
weapon  debris  Is  diffusion  controlled.  There  Is  also  evidence  that 
there  are  two  superimposed  distributions  of  soil  particles.  A 
computational  model  Is  developed  here  based  on  this  evidence.  The 
following  Is  a  synopsis  of  the  model:  A  portion  of  the  soil  Is  fully 
vaporized  along  with  the  weapon  debris  while  some  soil  Is  merely 
melted.  As  the  fireball  cools,  the  refractory  fission  products  will 
be  absorbed  by  this  liquid  material.  After  the  fireball  has  cooled 
below  the  soil  solidification  temperature,  the  remaining  fission 
products  can  be  adsorbed  onto  any  available  surfaces.  Any  soil  which 
entered  the  fireball  at  the  soil  solidification  time  or  later  will 
also  adsorb  fission  products.  Test  data  and  other  evidence  Indicate 
that  the  distributions  of  melted  and  unraelted  soil  particles  have 
different  modes.  This  model  uses  Henry's  Law  to  find  surface 
concentrations.  It  then  uses  diffusion  theory  to  transport  the 
fission  products  Into  the  particles.  In  addition.  It  allows  for 
Injection  of  unmelted  material  near  the  time  of  soil  solidification. 

The  results  of  the  research  Indicate  that  In  standard  DELFIC 
calculations  too  nwch  activity  Is  carried  In  the  larger  particles  and 


t.,  „„1.  tl-e  sMller  p.rt1cl«.  In  .ddUlon,  th.  distribution  o. 
«„.til.  fission  prod,«t  nuclides  rel.tlv.  to  .  refroctor,  reference 
^elioe  is  1n  generol  better  «de1ed  by  the  new  nethod.  Since  wny 
other  fellout  «del1n,  codes  «he  use  of  fits  to  OEtF.C  ectlrlty  size 
distributions,  these  codes  eight  be  .odifled  to  reflect  these  new 

findings. 


FALLOUT  FRACTIONATION  IN  SILICATE  SOILS 

I.  Introduction 

Background  r 

In  a  surface  burst  of  a  nuclear  weapon,  a  large  amount  of  soil  Is 
broken  up,  while  some  Is  melted  and  some  Is  even  vaporized.  Part  of 
this  debris  Is  taken  Into  the  fireball  and  Is  carried  with  It  as  the 
fireball  rises.  Part  of  this  material  will  come  Into  Intimate  contact 
with  the  radioactive  fission  products.  As  the  fireball  cools,  the 
soil  and  debris  from  the  device  will  condense  and  scavenge  the 
available  fission  products.  At  the  time  that  the  soil  begins  to 
solidify,  most  of  the  condensed  fission  products  will  be  more  or  less 
uniformly  distributed  throughout  the  particles.  After  that  time,  the  ' 
soil  particles  are  essentially  solid.  Those  fission  products  which 
were  still  gaseous  at  the  soil  solidification  temperature  will  later 
condense  onto  the  surfaces  of  the  particles  still  In  the  cloud. 

The  above  Is  a  very  simple  picture  of  the  processes  which  cause 
fractionation  taking  place  Inside  of  the  nuclear  cloud.  Simply 
stated,  radiochemical  fractionation  Is  a  distribution  of  Isotopes 
which  Is  different  from  that  which  would  be  expected  If  all  of  the 
radioactivity  were  distributed  uniformly  throughout  the  fallout 
particles.  An  Important  feature  of  this  process  which  has  always  been 
observed  but  neve-  treated  Is  the  presence  of  a  very  large  number  of 
Irregularly  shaped  particles  which  carry  radioactivity  only  on  their 


surfaces.  These  particles,  because  they  show  no  evidence  of  having 
been  oelted,  oust  have  been  Introduced  Into  the  active  region  of  the 
fireball  after  the  gas  temperature  had  dropped  to  near  or  below  the 
soil  solidification  temperature  (see  Figure  1).  One  of  the  objectives 
of  this  research  was  to  take  this  fact  Into  consideration. 

The  careful  treatment  of  this  problem  Is  very  Important  because 
of  the  large  difference  In  downwind  extent  of  Isodose  contours  If  the 
radioactivity  were  all  surface  distributed  as  opposed  to  being  all 
volume  distributed.  Figure  2  Illustrates  this  point  by  showing  In 
cross  section  the  dose  rates  along  the  hot  line  for  a  one  Mt  yield  (50 
percent  fission)  weapon.  In  this  case,  the  450  R  contour  for  the 
purely  volume  distributed  activity  extends  for  about  65  miles  past 
that  for  the  purely  surface  distributed  case.  This  Illustrates  that 
when  one  Is  trying  to  determine  the  expected  number  of  casualties  In  a 
given  attack  scenario,  the  effect  of  fractionation  cannot  be  Ignored. 
In  terms  of  placement  of  the  activity.  Figure  3  shows  the  distribution 
of  radioactivity  for  these  two  extremes  as  a  function  of  particle  size 
for  an  assumed  lognormal  particle  size  distribution. 

Historically,  the  first  attempt  to  develop  a  theory  for 
fractionation  was  by  Miller  with  his  Thermodynamic  Equilibrium  Model 
(Miller,  1960;  Miller,  1963;  Miller,  1964).  The  model  essentially 
consisted  of  distributing  the  fission  products  In  the  cloud  among  the 
particles  according  to  the  equilibrium  distribution.  The  model  used 
1400  degrees  C  as  the  temperature  below  which  condensing  nuclides 
could  not  penetrate  the  fallout  particles.  The  difficulty  with  the 
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Figure  1.  Schematic  of  Early  Cloud  Near  the  Time  of 
Soil  Solidification 
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Figure  2.  Surface  versus  Volume  Distributed  Activity 
for  a  1  Megaton  Ground  Burst 
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■odsi  was  that  adequate  thermodynamic  data  was  not  available  at  the 
time  It  was  developed.  Thus,  the  model  could  not  be  properly 
1nq)lemented.  The  next  attempt  to  treat  fractionation  came  from 
Frelling  with  h1$  Radial -power-distribution  Model  (Frelllng,  1961; 
Frelllng,  1963;  Frelllng,  1963  b;  Frelllng,  1964).  This  model  assumed 
that  there  was  a  collection  of  spherical  particles  and  that  all  mass 
chains  were  distributed  according  to  some  power  of  the  particle 
radius.  With  these  assumptions,  Frelllng  found  that  the  model 
produces  logarithmic  correlations  between  radionuclide  ratios. 
Frelllng,  In  his  analysis  of  Pacific  weapons  test  data,  observed  that 
radionuclide  ratios  could  be  correlated  logarithmically  as  well  as 
they  could  be  linearly.  This  was  taken  as  evidence  that  the  method 
had  some  validity.  The  correlation  parameters  formed  the  basis  for  a 
model  to  predict  fractionation  effects.  The  next  attempt  to  treat  the 
problem  was  by  Korts  and  Norman  with  their  diffusion  limited  approach 
(Korts,  1967).  In  developing  this  model,  a  great  deal  of  work  was 
done  to  measure  the  thermodynamic  properties  of  fission  product  oxides 
and  the  diffusion  properties  of  various  soil  types.  The  essence  of 
the  method  was  a  hybrid  of  the  equilibrium  treatment  of  Miller  with 
diffusion  of  the  condensed  fission  product  oxides  Into  the  soil 
particles.  This  model  assumed  a  distribution  of  spherical,  glassy 
particles.  Attempts  were  underway  In  the  late  1960 's  to  bring  even 
greater  sophistication  to  the  fractionation  problem,  the  skeleton  for 
a  more  complete  kinetic  model  had  been  laid  out  by  this  time,  but 
because  of  funding  cutbacks  and  laboratory  reorganizations  (e.g.,  the 
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closing  of  the  U.S.  Naval  Radiological  Defense  Laboratory),  research 
on  particle  formation  was  largely  stalled. 

Problem  Statement 

The  objective  of  this  research  project  was  to  review  the 
literature  on  fallout  fractionation  and,  if  possible,  extend  or 
Improve  the  theory.  As  a  result  of  ^e  research,  a  complete 
computational  model  was  produced  even  though  that  was  not  one  of  the 
initial  goals.  In  fact,  two  distinct  models  were  produced:  one  is 
referred  to  as  the  G-X  uiodel  (for  glass-crystal),  while  the  other  (a 
hybrid  of  Miller's  model)  is  referred  to  hereafter  as  the  S-Y  model 
(for  surface-volume).  The  latter  model  was  incidental  to  the  main 
thrust  of  the  research.  For  this  reason  it  is  presented  in  Appendix  A 

Scope 

The  problem  of  fractionation  includes  many  aspects  which  will  not 
be  treated  here.  Mo  attempt  will  be  made  to  treat  fractionation  which 
occurs  because  of  sampling  methods.  Nor  will  fractionation  which 
occurs  because  of  leaching  of  grounded  particles  prior  to  collection 
be  treated.  The  treatment  here  is  of  the  process  termed  "primary 
fractionation"  by  Freiling  (Freiling,  1961:  1991).  The  problem  is 
further  limited  to  surface  bursts  over  silicate  soil.  This 
restriction  to  silicate  soils  is  made  since  no  effort  has  ever  been 
made  to  accurately  identify  the  particle  sire  distributions  for  the 
different  types  of  particles  found  in  surface  bursts  over  calcium 
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based  sofis.  The  restriction  to  surface  bursts  Is  made  because  the 
fractionation  process  Is  fundamentally  different  for  each  type  of 
burst  condition.  Surface  bursts  are  of  particular  Interest  because 
they  present  the  greatest  fallout  hazards.  They  are  also  Interesting 
because  of  the  complicated  mixing  of  unmelted  material  with  fission 
products,  melted  soil  and  weapon  debris. 
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II.  HIstoHcjal  Development 


Introduction 

The  early  attempts  to  model  fractionation  were  made  difficult  by 
a  general  lack  of  data.  The  collection  of  v^eapons  test  data  which 
would  be  of  use  to  fallout  researchers  was  a  low  priority.  In 
addition,  early  Investigators  were  not  sure  what  types  of  data  were 
needed.  Oa:a  collection  for  fallout  studies  was  often  an 
afterthought.  Even  when  it  wa$  carefully  planned,  equipment  was  often 
unreliable.  High  radiation  in  some  areas  would  frequently  delay 
sample  recovery  allowing  contamination  of  samples  by  non- radioactive 
debris  carried  by  winds  {Freiling,  1965  b). 

On  top  of  these  difficulties,  nuclear,  thermodynamic,  and 
chemical  data  for  the  fission  products  and  their  oxides  were  limited 
or  non-existent.  In  spite  of  these  problems,  several  investigators 
managed  to  develop  workable  models  for  predicting  the  effects  of 
fractionation  as  well  as  a  theoretical  framework  for  understanding . 
it.  The  most  notable, of  these  were  E.  C.  Freiling,  J.  H.  iJorman,  and 
C.  F.  Miller.  In  the  following  sections,  their  models  wi ,  be 
surveyed  in  order  to  provide  a  foundation  for  understanding  the  6-X 
model  developed  in  the  next  chapter. 

Miller's  Thermodynamic  Model 

Introduction.  Miller's  theory  of  fractionation  was  developed 
about  the  same  time  as  Frei ling's  treatment  (Miller,  1963).  Since 


Miller's  model  required  rather  a  great  deal  to  be  known  about  the 
thermodynamics  of  the  cooling  fireball,  and  since  the  thermodynamic 
data  for  the  fission  products  v/ere  not  available,  it  was  of  limited 
value  as  a  computational  model.  It  did  in  its  original  form,  however, 
treat  many  of  the  important  phenomena  in  a  qualitative  way  and  pointed 
the  way  for  the  collection  of  the  data  that  would  be  necsssar:'  to 
implement  a  more  complete  theoretical  model  at  some  later  date.  It 
was  because  of  the  lack  of  data  for  Miller's  model  that  Freiling's 
semi -empirical  Radial  Power  Law  was  adopted  as  a  stop-gap  device  for 
predicting  fractionation  effects  (Freiling,  1965).  But  even  as  a 
purely  theoretical  model,  Miller's  treatment  leaves  room  for 
improvement  as  it  neglects  some  Important  phenomena. 

It  should  be  noted  in  the -foil owing  discussion  that  Hiller  treats 
the  melt  of  carrier  material  as  a  single  mass  rather  than  as  a 
distribution  of  particles  of  varving  sizes.  This  poses  no  particular 
problem,  although  it  Ignores  the  effect  of  particle  size  on  the 
condensation  mathematics. 

Miller  divides  the  condensation  process  into  tv/o  more  or  less 
distinct  time  periods.  The  major  feature  of  the  first  period  is  the 
existence  of  vapor-liquid  phase  equilibria.  This  period  of 
condensation  ends  when  the  carrier  material  solidifies,  with  the 
fission  products  either  fixed  in  a  solid  solution  matrix  or  compounded 
with  the  carrier  material. 

The  major  feature  of  the  second  period  of  condensation  is  the 
existence  of  vapor-solid  phase  equilibria  in  which  the  remaining 
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fission  product  elements  condense  at  lower  temperatures  on  the 
surfaces  of  solid  particles.  The  second  period  of  condensation  never 
ends  completely  except  for  those  particles  which  leave  the  volume  of 
space  containing  the  residual  gases.  In  fact>  the  process  can  reverse 
for  a  fission  product  element  that  later  decays  to  a  more  volatile 
element.  For  example,  elements  like  Iodine  and  the  rare  gases  could 
sublime  as  fast  as  they  form  from  non-volatile  precursors  which 
condensed  on  the  surfaces  of  fallout  particles.  This  process  Is 
unlikely,  however,  when  the  fission  products  are  trapped  within  a 
glassy  iwtrlx.  The  vapor  pressures  due  to  the  low  concentration  of 
dissolved  fission  products  would  be  extremely  low,  and  diffusion 
through  the  solid  glass  would  be  very  slow. 

The  essential  problem  In  the  theory  for  the  process  during  the 
first  period  of  condensation  Is  to  establish  the  vapor-liquid  phase 
equilibria  of  each  fission  product  element  at  the  time  that  the 
carrier  material  solidifies;  that  Is  to  determine  the  fraction  of  each 
element  present  which  has  condensed  and  migrated!  Into  the  carrier  melt 
at  the  time  of  solidification. 

When  one  of  the  two  phases  In  contact  Is  a  gas,  simple  kinetic 
theory  can  be  used  to  r; how  that  condensation  -  vaporization 
equilibrium  can  be  established  very  quickly  at  temperatures  above  2000 
degrees  Kelvin  (Miller,  1960).  Thus,  those  gaseous  species  of  each 
fission  product  element  that  do  not  react  with  the  liquid  carrier  but 
dissolve  into  It  should  obey  Henry's  Law  of  dilute  solutions  (see 
Figure  4).  In  fact,  the  solutions  should  be  sufficiently’  dilute  as  to 
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Ffgure  4.  Henry’s  Law  of  Dilute  Solutions 
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result  in  no  appreciable  change  in  the  Gibbs  free  energy  of  the  liquid 
carrier.  Thus,  the  free  energy  of  each  element  in  the  solution  should 
be  Independent  of  any  other.  It  is  possible,  then,  to  consider  the 
solubility  of  each  element  by  treating  each  as  a  binary  system  with 
the  carrier.  Further,  there  should  be  no  appreciable  surface  loading 
(large  excess  surface  concentrations  of  fission  product 
elements/compounds)  during  the  condensation  process  if  the  temperature 
range  over  which  the  liquid  carrier  exists  exceeds  200  or  300  degrees 
C.  Miller  allows,  however,  that  a  concentration  gradient  should  exist 
for  the  larger  particles  of  which  some  may  not  be  melted  in  the  center 
when  the  air  or  gas  temperature  about  the  particle  falls  below  the 
melting  point  of  the  bulk  carrier  (Miller,  1963:  99>101). 

Development  of  Equations  for  Hiller’s  Model.  Hhen  considering 
the  thennodynamlcs  of  a  multicomponent  system,  the  various 
thermodynamic  functions  must  include  the  amount  of  each  component  as  a 
variable.  The  chemical  potential,  u^,  is  a  useful  quantity  for 
treating  equilibrium  or  near  equilibrium  systems.  It  represents  the 
change. In  the  Gibbs  free  energy  of  a  system  per  mole  of  added 
component  i,  with  te^erature,  pressure,  and  the  other  molar 
quantities  kept  constant.  Adamson  gives  the  total  derivative  of  the 
chemical  potential  as  (Adamson,  1373:  355): 

d«^  -  RT  d(lnf^)  (1) 
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where  Is  the  fugacity  of  component  i.  Fugacity  Is  an  effective 
pressure  fcr  the  component.  It  accounts  for  deviations  from  Ideal 
behavior.  The  defining  equation  for  fugacity  for  state  1  «  1  going  to 
state  1  »  2  Is  (Andrews,  1970:  204): 


,  A  ®2  .  11m 


Pl>0  '’l 


W 


(2) 


where  Is  the  Gibbs  free  energy  for  state  1.  When  dX^  moles  of 
element  1  condense  from  a  gaseous  mixture  to  a  liquid  solution  leaving 
moles  In  the  qas  and  X^  moles  In  the  liquid,  the  change  In  the 
chemical  potential  In  the  gas  Is  found  as  follows:  First  we  expand 
the  derivative  given  In  Equation  (1)  to  obtain 


•  ^^1  - 

d.i .  ‘^P,x;  ^ 


(3) 


or  equivalently 


.  »lnf.  alnf.  »lnf. 

<■“)  ■  ‘  "T'-jTJ-jdT  . 

ifX^  P,X^  3X^  T,P  ^ 


(4) 


The  corresponding  change  In  the  chemical  potential,  U|,  In  the 
solution  Is 
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(8) 


•nd  Equation  (5)  becomes 


dii^  •  dP  “  3^1  dT  *  RT 


dXi 


(13) 


The  tMO  changes  In  the  chemical  potential  for  the  transfer  of 
dX^  moles  from  the  gas  mixture  to  the  liquid  solution  are  equal  In 
an  equilibrium  process.  Dividing  Equations  (12)  and  (13)  by  RT  while 
noting  that 


(14) 


and  then  equating  Equation  (1)  and  Equation  (2),  after  substituting 
for  the  Indicated  partial  differentials,  gives 


^  ’ll  .  ’T) 

KT  <IP  -  — 2  if  ♦  d  Inx  .  ^  (JP  -  — J  dT  *  d  In  X,  (IS) 
RX 


RT' 


In  which TJ*  Is  the  partial  molar  volume  of  element  1  In  the  gas 
mixture,  ^  Is  Its  relative  partial  molar  heat  content  In  the  gas 
mixture,  Is  Its  partial  molar  volume  In  the  liquid  solution. 

Is  Its  relative  partial  molar  heat  content  In  the  liquid 
solution,  R  Is  the  molar  gas  constant,  P  Is  total  pressure,  .'nd 
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T  1$  the  temperature  In  degrees  Kelvin.  In  the  case  of  uTiuxe 
solutions  (Miner,  1960:  10) 

TjJ  -  .  A  (16) 

where  aH^  Is  the  heat  of  vaporlaatlor.  of  the  condensing  element  1. 
Miller  points  out  that  aH^  Is  the  heat  of  reaction  for  the 
vaporization  of  the  gas  from  solution.  If  It  exists  as  a  different 
compound  In  solution,  the  heat  of  fonnatlon  of  this  compound  Is 
Included  In  the  value  of  aH^.  For  an  Ideal  solution,  aH^  Is  just 
the  heat  of  vaporization.  Since 

^1 

7^  »  T  then  {^)  (17) 

can  be  neglected.  Thus,  Equation  (15)  becomes 
V*  aH 

^  dP  ♦  ddnX!)  .  — 7  dT  ♦  ddn  X^) 

-  (18) 

but  by  Equation  (10)  this  reduces  to 

aH 

^  ♦  ddnXp  »  — 7  dT  ♦  ddn  X^) 

F  ’  RT  (191 

Upon  Integration,  this  becomes 
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In  P  *  In  '  *  In  *  In 


(20) 


or 


and  finally 


AH^ 

KT 


(21) 


P  x!  -aH^/rt 
Tf  •  •'i  * 


(22) 


o 


In  which  Is  an  Integration  constant  and  where  the  term 


k,  S  k;  «  * 


(23) 


can  be  Identified  from  Equation  (€)  and  Equation  (8)  where  (f^  « 
f^  »  p^,  the  partial  pressure  of  element  1;  and  f^  ■  P)  as  the 
Henry's  Lav/  constant.  Now  substituting  from  the  definitions  of 


and 


X*  *  "l  "l 

^  "I  *  "air  "air 


(24) 
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(25) 


then  Equation  (22)  becomes 


"l/"a1r  •'l 

"(^"soir  ^ 


(27) 


where  n^^^^  Is  the  number  of  moles  of  liquid  soil.  Finally  (from 
Equation  27),  the  number  of  moles  of  the  1th  element,  n^.  In  the 
liquid  Is  related  to  the  total  number  of  moles  In  both  phases  by 
the  relation 


(28) 


In  the  second  period  of  condensation  under  the  Miller  model,  the 
soil  carrier  has  solidified.  Fission  product  uptake  then  continues  by 
surface  adsorption.  The  adsorption  condensation  can  be  considered  as 
one  In  which  the  relative  amount  of  each  element  condensed  Is  related 
to  Its  sublimation  pressure  (Miller,  1960:  109-110).  The 
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computattonal  values  from  such  a  process  will  reflect  the  relative 
volatility  of  the  constituent  molecules  at  all  temperatures  at  which 
this  kind  of  condensation  can  occur.  If  an  excess  of  solid  surface 
area  is  present,  the  number  of  moles  condensed  by  the  process 
(assuming  that  the  process  of  reversible)  at  ar\y  time  after 
solidification  of  the  carrier  Is  given  by 


I 


(29) 


In  which  n^  Is  the  residual  amount  of  element  1  condensed  on  the 
surface  of  the  solid  particles,  n^  Is  the  amount  of  the  element  In 
the  gas  phase  at  the  start  of  the  second  period,  and  n!^  Is  the 
residual  amount  In  the  vapor  phase.  All  three  quantities  depend  on 
time  because  the  radionuclides  continue  to  decay.  In  the  case  of  a 
perfect  gas 


P*v 
"l  •  “ITT 


(30) 


where  v  Is  the  volume  enclosing  the  n^  moles  of  the  gaseous  species 
and  Is  the  sublimation  pressure  and  Is  given  by 


-  aS./R  -aH  /RT 
Pj  -  e  *  e  * 


(31) 


.  I 

1 

I 

I 

I 
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In  which  Is  the  entropy  of  sublimation  and  Is  the  neac  of 
sublimation  at  the  temperature,  T.  If  mixing  In  the  fireball  volume 
Is  assumed  to  be  uniform  for  the  particles  and  gaseous  fission 
products,  then  v  Is  the  fireball  volume.  The  material  balance  for 
ele.nent  1  Is 


•  •  "I 


(32) 


so  that.  Equation  (29),  with  Equations  (27),  (28),  (30),  and  (32), 
becomes 


In  Equations  (28)  and  (33),  the  quantities  n^Q^-j  and  v  ere 
estimated  from  the  thermal  data  and  empirical  scaling  equations. 

Miller's  work  was  notable  In  that  It  treats  the  essence  of  the 
problem  1n  a  physically  meaningful  way.  The  method  developed  by 
Frelling  presented  In  the  next  section  Is  largely  emplHcal.  While 
none  of  the  necessary  Henry's  Law  constants  were  available  when  Miller 
wrote  his  report,  many  are  now.  In  those  cases  where  data  Is  still 
not  available,  Raoult's  Law,  which  Is  an  Idealization  of  Henry's  Law, 
may  be  used.  In  Appendix  A,  a  variation  of  the  Miller  method  Is 
developed  which  has  an  especially  attractive  feature:  viz..  It  could 
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be  easily  modified  for  use  with  fast  running  "smear  codes"  such  as 
that  developed  by  Bridgman  and  Bigelow  (Bridgman,  1982). 

The  next  section  provides  an  overview  of  the  Radial  Power  Law 
developed  by  FrelHng.  While  It  Is  not  as  elegant  theoretically  as 
Miller's  equilibrium  treatment.  It  has  the  advantage  of  being  less 
difficult  computationally. 

RADIAL  POWER  LAW 


Introduction.  FrelHng's  Radial  Power  Law  Is  a  semi -empirical 
fractionation  model  which  Is  presently  used  In  modified  form  In  the 
official  Department  of  Defense  fallout  modeling  computer  code, 

DELFIC.  The  discussion  here  Is  restricted  to  the  pure  form  of  the 
Radial  Power  Law  rather  than  the  modification  used  In  DELFIC. 
Discussion  of  the  Tompkins  Implementation  of  the  Radial  Power  Law  In 
DELFIC  form  Is  presented  In  the  following  section.  The  Radial  Power 
Law  has  some  limitations,  which  will  be  discussed  after  the  model  has 
been  presented. 

The  model  Itself  Is  simple  and  does  to  a  degree  Incorporate  the 
radiochemistry  of  the  cooling  radioactive  elements.  First,  the 
principal  compound  (usually  but  not  always  an  oxide)  of  each  element 
which  would  exist  In  the  hot  fireball  must  be  Identified.  The  boiling 
point  of  this  compound  Is  tht'n  taken  as  the  boiling  point  of  that 
element  for  computational  purposes.  Next,  a  melting  point  is 
determined  for  the  soil  which  will  be  the  carrier  of  the 


rad^oact^vity.  Without  the  presence  of  a  large  mass  of  this  carrier 
material,  there  Is  essentially  no  local  fallout  from  a  nuclear  weapon 
(Miller,  1963:  3).  A  rough  rule  of  thumb  Is  that  particles  larger 
than  10-20  urn  will  be  deposited  locally,  or  within  a  few  hundred  miles 
of  the  detonation  point  (Glasstone,  1967:  36-371.  In  the  Freillng 
model,  the  user  must  supply  the  yield,  height  of  burst,  soil  type, 
fuel,  fission  yield,  the  particle  size  distribution,  and  the  fission 
yields  of  each  radionuclide  for  each  of  several  fuel  types  and  neutron 
energy  spectra.  In  additicn,  the  model  needs  all  of  the  decay  chain 
data  necessary  to  follow  the  decay  of  the  fission  products. 
Thermodynamic  characteristics  of  the  fireball  must  be  computed  as 
well.  These  Include  how  much  soil  .  vaporized  or  melted  and  carried 
aloft,  as  well  as  the  time  at  which  the  fireball  reaches  the 
solidification  temperature  of  the  soil.  The  number  of  fission  product 
atoms  which  are  In  a  refractory  form  at  that  point  In  time  Is 
determined  for  each  mass  chain  by  solving  rigorously  the  Bateman 
equations  which  describe  radioactive  decay.  The  fraction  of  atoms  In 
each  chain  which  are  In  a  refractory  form  at  this  time  Is  used  as  a 
parameter  In  distributing  the  nuclides  of  each  chain  In  the  soil 
particles. 

D1 scusslon.  There  are  two  essential  points  In  the  following 
discussion:  F^rst,  Frelling's  conclusions  from  his  analysis  of 
Pacific  weapons  test  data  are  not  well  supported  by  data  from  ground 
bursts  over  silicate  soil.  Second,  Frelling's  theoretica-  explanation 


of  Ms  empirical  correlations  Is  likewise  not  well  supported  by  later 
data.  The  data  from  continental  United  States  tests  suggest  that 
simple  empirical  correlations  are  not  sufficient  to  predict 
fractionation  behavior  In  the  case  of  ground  bursts  over  silicate 
soils. 

Frelling  makes  the  following  analysis:  In  a  sample  of  fallout, 
the  apparent  number  of  fissions  which  yielded  a  given  nuclide  Is 

f^  -  a^/Y^  (34) 

where  a^  Is  the  number  of  atoms  of  nuclide  1  In  the  sample  and 
Is  the  number  of  atoms  of  nuclide  1  per  fission.  The  ratio  r^  yis 
defined  by 

This  ratio  Is  useful  for  the  following  reason:  Let  nuclide  j  be  a 
nuclide  which  Is  volatile  and  whose  precursors  were  also  volatile. 

This  nuclide  would  tend  to  be  surface  distributed  since  both  It  and 
Its  parents  would  have  lower  boiling  points  than  the  soil  carrier. 

Now  let  nuclide  1  first  be  another  similarly  behaving  nuclide.  Then 
r^  j  would  be  approximately  equal  to  1.  Now,  If  j  was  a  nuclide 
which  was  refractory  and  had  refractory  parents,  then  the 
corresponding  boiling  points  would  be  higher  than  the  soil 
solidification  temperature.  So  the  tendency  would  be  for  this  nuclide 
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to  be  volume  distributed.  The  ratio  j  would  then  have  an 

extremum  since  we  are  measuring  the  greatest  possible  difference  In 

behavior  of  two  nuclides.  How  we  take  r  ,  to  he  our 

''ref  *ref 

reference  r^  j  where  v^^  Is  the  reference  volume  distributed, 
refractory  nuclide  and  s^^  Is  the  reference  surface  distributed, 
volatile  nuclide.  This  ratio  will  be  abbreviated  r„  .  (an  alternate 

¥»S 

notation  Is  r  where  the  subscript  r  Indicates  a  refractorily 
behaving  nuclide  and  v  represents  a  volatllely  behaving  nuclide). 

Then  form  the  ratio  r^  ^  between  f^  where  1  Is  arbitrary  and  f^ 
where  s  represents  our  reference  surface  distributed  nuclide.  Next, 
consider  plotting  the  relationship  between  these  ratios  for  various 
actual  fallout  samples.  That  Is,  given  a  sample  of  fallout  material, 
analyze  the  characteristic  end  products  of  each  mass  chain  (®®Sr  for 
1  ^  90,  for  1  .  99,  ^^^Cs  for  1  -  137,  etc.).  Then  form  the 
ratios  r^  ^  and  r^  ^  (we  will  take  s  ■  89  for  ®®Sr  and  v  ■  95 
for  55zr  for  our  reference  nuclides).  These  ratios  can  best  be 
correlated  logarithmically  for  the  high  yield  coral  and  water  surface 
bursts  analyzed  by  Frelling  (Frelling,  1961:  1994,  Heft,  1970: 
256-257).  Frelling  considered  two  types  of  relationships  between 
fractionation  ratios: 

’^1,89  ■  *^1  ’'95,89  *  *^2 

(b)  In(r^^gg)  «  kj  "^nlr^g  ♦  k^ 
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(36) 

(37j  . 


S  *•! 


^  (') 


where  the  k's  'ere  constants  to  be  determined.  Of  these  two,  the 
latter  had  the  smaller  variance  In  the  data  fit  (Note:  «nd 

Tjg  gg  are  functions  of  particle  size).  Note  further  the 
following:  If  1  Is  a  refractorily  behaving  nuclide,  then  r^^gg  will 
mimic  rgg  gg  so  that  Equation  (37)  reduces  to 

kg  >  (ln(ry  j)  -  k^)/ln(ry^j)  (38) 


If  k^  was  zero,  then  we  see  that  kg  ■  1.  Similarly,  If  1  Is 
volatile,  then  Equation  (37)  reduces  to 


In  j  .  kj  In  r,  ^  ♦  k^ 


but  r.  .  ■  1  so  we  have 

S  f  9 


0  .  kg  In  r^^j  ♦  k^ 


*^3  • 


and  again  If  k^  «  0,  ihen  kg  «  0.  As  It  turns  out,  Frelling  found 
In  his  plots  that  k^  Is  approximately  zero  and  kg  Is  approximately 
unity  when  1  Is  a  nuclide  whose  chain  Is  all  refractory,  and  kg  Is 
approximately  zero  when  1  Is  a  nuclide  whose  chain  Is  all  volatile. 
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In  sunniary,  logarithmic  correlations  seem  to  work  we11  for  high  yield 
coral  surface  bursts.  They  do  not,*  however,  for  low  yield  silicate 
surface  bursts.  Later  Investigators  found  evidence  for  linear 
relationships  In  the  latter  case  (Heft,  1970:  256,  Bridgman,  1983). 

To  establish  a  theoretical  basis  for  his  semi -empirical  model, 
Frelling  postulated  first  that  the  particle  size  distribution  was 
lognormal.  He  then  further  postulated  that  the  distribution  of  a 
given  mass  chain  goes  as  some  power  of  the  radius  of  the  fallout 
particles,  where  the  exponent  lies  between  2  and  3  (Tompkins,  1968: 
11-16).  If  the  chain  Is  all  refractory,  then  Its  long-lived  members 
will  be  refractory  and  will  be  distributed  as  the  third  moment  of  the 
particle  size  distribution.  But  If  the  chain  Is  all  volatile,  then 
the  members  observed  1n  analysis  will  be  volatile  and  will  bi:  , 
distributed  as  the  second  moment  of  the  final  particle  size 
distribution.  When  Frelling  used  this  type  of  analysis,  he  found  that 
the  data  for  coral  and  water  high  yield  surface  bursts  best  fit 
Equation  (37)  with  k^  «  0  and  kj  equal  to  the  square  root  of  the 
fraction  of  the  decay  chain  1  which  exists  In  refractory  form  at  the 
time  of  soil  solidification  (Frelling,  1961:  1995).  This  conclusion 
violates  the  Intuitive  hypothesis  that  volatile  nuclides  are  surface 
distributed  and  refractory  nuclides  are  volume  distributed.  If  It 
were  true  that  a  single  mass  chain  contained  some  nuclides  which  were 
purely  volume  distributed  and  some  nuclides  which  were  purely  surface 
distributed,  then  the  superposition  of  these  two  distributions  would 
not  be  a  simple  fractional  moment  between  2  and  3  of  the  lognormal 
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particle  size  distribution  (which  was  Itself  a  lognormal 
distribution).  It  rather  Is  the  sum  of  two  different  lognormal 
distributions,  which  Is  not  lognormal. 

To  help  resolve  the  apparent  conflict  between  Frel ling’s  theory 
and  the  data,  consider  the  following  points:  First,  the  data  used  by 
Frelling  was  limited  to  high-yield  surface  oursts  over  coral  and 
water.  There  Is  evidence  that  for  silicate  soil  surface  bursts  of 
lower  yield  the  slope  should  be  equal  to  the  fraction  of  the  chain 
which  exists  In  refractory  form  at  the  time  of  soil  solidification 
rather  than  the  square  root  of  the  refractory  fractlon.vT^,  as 
Fre111ng*s  Pacific  Test  Site  data  suggested  (Bridgman,  1983:  b)  (see 
Figures  S  and  6).  So  Frelling's  square-root  relationship  may  only 
apply  In  an  approximate  sense  to  high  yield  coral  or  water  surface 
bursts. 

Second,  when  the  soil  carrier  Is  In  the  liquid  phase,  the  degree 
to  which  the  radioactive  nuclides  collect  on  and  diffuse  Into  the  soil 
droplets  Is  a  function  of  the  vapor  pressure  of  those  nuclides.  As 
noted  1n  the  previous, section.  Miller  asserts  that  the  essence  of  the 
fractionation  problem  lies  In  vapor-liquid  phase  equilibria 
established  for  each  fission  product  element  (or  compound)  at  the  time 
Just  prior  to  solidification  of  the  carrier  material.  That  Is,  one 
needs  to  find  what  fraction  of  each  fission  product  could  be  found  In 
the  liquid  phase  and  determine  to  what  degree  It  had  diffused  Into  the 
carrier  melt  when  solidification  occurred.  Those  gaseous  species  of 
each  fission  product  element  that  do  not  react  with  the  liquid  carrier 
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OBSERVED  SLOPE 


Figure  5.  The  Square  Root  of  versus  the  Correlation  Slopes 

for  Johnnie  Boy  Data 


OBSERVED  SLOPE 


Figure  6.  versus  the  Correlation  Slopes 
for  Johnnie  Boy  Data 


3 


but  dissolve  Into  It  should  obey  Henry's  Law  of  dilute  solutions.  A 
concentration  gradient,  however,  should  exist  In  many  particles, 
especially  for  the  larger  particles  of  which  some  not  have  been 
completely  melted. 

In  other  words,  one  should  not  expect  to  find  all  refractory 
nuclides  distributed  unlfomnly  throughout  the  volume  of  the  soil 
particles.  Similarly,  volatile  nuclides  will  not  necessarily  be 
restricted  to  the  surfaces  of  the  soil  particles.  On  the  other  hand, 
there  exists  the  possibility  for  desorption  of  nuclides  which  have 
condensed  on  the  surfaces  of  particles,  especially  after  they  have 
left  the  cloud.  Some  nuclides  will  have  especially  small  sticking 
coefficients  and  thus  exhibit  poor  surface  adsorption.  The  lack  of 
treatment  ot  these  phenomena  Is  a  weakness  of  the  Radial, Pow^r  Law 
model . 

Mow  consider  the  following  analysis.  First,  define 

d  >  particle  diameter 

Mid)  a  particle  size  distribution 

Fj  a  total  number  of  equilavent  fissions  In  all  particles 

,  a  fission  yield  of  1-th  mass  class 

F^{d)  a  distribution  of  equivalent  fissions  for  particles 
of  diameter  d  for  mass  chain  1 
(r,  ,)  a  fractionation  ratio  n  F<(d)/F.(d)  where  r  Is 
Index  for  a  refractory  behaving  chain 
fj.  a  fraction  of  the  decay  chain  1  which  exists  in 

refractory  form  at  the  time  of  soil  solidification 
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While  FreHIng  used  a  volatile  chain  for  his  reference  chain,  a 
refractory  chain  Is  the  more  natural  choice  since  It  will  be  more  or 
less  uniformly  distributed  throughout  the  fallout.  The  following 
analysis  will  therefore  use  a  refractory  chain  for  the  reference 
nuclide.  Now  If  1t*s  true  that  the  refractory  nuclides  will  be 
distributed  as  the  third  moment  of  the  particle  size  distribution  and 
volatile  nuclides  will  be  distributed  as  the  second  moment,  then  for  an 
arbitrary  chain  1  we  can  write 

"^N(d)  |fpd^  ♦  {l-f^)d^[d  d  (42) 

where  and  are  normalization  constants.  For  a  refractorily 
behaving  chain  we  get 

F,  N(d)d^ 

F  (d)  .  , 

(43) 

Now  the  fractionation  ratio  Is 


)(d) 


Fyf?(d)  (f^d^  ♦  (l-%)d^) 
{N(d)  (f^d^  ♦  {l-f^)d^  )d  d 
F^  N{d)d^ 
jjM{d)d^  d  d 


(44) 
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Now  notice  In  Equation  (48)  the  expression 


■  • 

■  f  i 


/yd){f^  ♦  (l-f^)/d  [d  d 


(49) 


Is  a  normalization  constant  for  mass  chain  1.  For  a  chain  which 
consists  only  of  volatile  species.  Equation  (48)  becomes  then 


(r^^^)(d)  - 


since  f^  •  0  In  this  case.  While  for  a  general  chain 
(ri  r)(d)  -  C^|f^  ♦  (l-V/d[ 


(50) 


(51) 


Solving  Equation  (50)  fof  d  and  substituting  Into  Equation  (51)  yields 


(52) 


This  Is  a  simple  linear  relation  between  (r^  .)(d)  and  (r„  ,)(d). 

Now  Frelling  suggests  that  Instead  of  separating  the  surface  and 
volume  fractions,  that  a  given  mass  chain  will  have  Its  activity 
distributed  as  a  simple  power  of  the  diameter  of  the  particles.  The 
equations  equivalent  to  Equation  (50)  and  Equation  (51)  under 
Frel ling's  assumption  (the  Radial  Power  Law)  are 
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(53) 


(r,_r)(d) 


and 


C|  d 


(54) 


Eliminating  d  between  these  last  two,  and  talcing  the  logarithm  of  both 
sides  of  the  result,  gives 


ln(r^^^)(d)  -  (l-v?j;)  ln(r^^^)(d) 
♦  i(v77-l)ln  C  In  C 


(55) 


This  Is  a  logarithmic  relation  with  slope  (1  -  vTp,  Frelllny  has 
observed  this  logarithmic  relationship  In  fallout  sample  data,  and 
Indeed  the  1  -  expression  In  the  slope  above  Is  a  best  fit  of 
the  observed  correlation  slopes  to  the  fraction  of  the  chain  existing 
In  refractory  form  at  the  time  of  soil  solidification  (Frelllng,  1961: 
11-16). 

Note  that  If  Equation  (52)  Is  plotted  on  log-log  paper  the  limits 
(fj.  «  0,  fj.  »  1)  agree  with  observations  In  the  data,  but  for 
Intermediate  chains  (0  <  f ^  <  D.  Equation  (52)  predicts  a  distinct 
curvature  to  the  data  (see  Figure  7).  This  doesn't  seem  to  agree  very 
well  with  the  data  cited  by  Frelllng.  On  the  other  hand,  the 
expression  In  Equation  (55)  agrees  very  well  with  the  Pacific  weapons 
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Fi'gure  7.  Plots  of  Equation  (52)  for  Several  Values 
of  the  Frelllng  Ratio 
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test  data.  Russell  and  Heft  (Heft,  1970:  2E6-257,  Bridgman,  1983  c) 
have  shown,  however,  that  1n  the  case  of  silicate  soil  surface  bursts, 
that  we  do  see  curvature  in  the  plots  of  Equation  (52). 

It  is  clear,  then,  that  there  are  limitations  to  the  Freiling 
model.  Freiling  himself  offered  several  caveats  on  the  method 
(Freiling,  1961:  1958,  Radionuclides,  1970:  336).  The  correlations 
that  Freiling  has  observed  should  depend  on  weapon  yield,  and  the  type 
of  carrier  which  constitutes  most  of  the  particles.  Yet  his 
corre’ations  involved  only  high  yield  data  for  coral  and  water  surface 
bursts.  Since  Freiling  made  his  correlations  irrespective  of  yield  or 
carrier  material,  the  Radial  Power  Law  disregards  these  effects.  His 
approach  is  primarily  a  graphical  device  which  is  useful  for  drawing 
general  trends  in  behavior.  A  single  set  of  correlations  cannot  be 
expected  to  cover  all  burst  conditions.  Care  should  be  exercised  in 
applying  the  Radial  Power  Law  to  low  yield  events  or  to  events  over 
silicate  soils. 

In  the  next  section,  the  Tompkins  modification  to  the  Radial 
Power  Law  will  be  discussed.  It  was  developed  to  bring  the  results  of 
the  Radial  Power  Law  into  better  agreement  with  specific  activity 
observations. 

MODIFIED  RACIAL  POVIER  LAW 

Introduction.  When  the  Radial  Power  Law  was  considered  by 
Defense  Atomic  Support  Agency  (DASA)  for  the  DELFIC  code,  Tompkins 


*-  JT;  *t“rf  IT*  • -« 


V^^iiTCT iTJ.TWT'T.’  '•r.’  yr^  if ^  ir^ 


made  the  observation  that  there  is  a  tendency  for  specific  activity 
data  to  level  out  as  particle  size  Increases  regardless  of  the 
volatility  of  the  species  being  Investigated.  The  transition  point  In 
the  curves  appeared  to  Tompkins  to  be  about  100  to  200  microns  for 
surface  bursts  (Tompkins,  1968:  15).  In  order  to  account  for  this, 
Tompkins  required  all  fission  products,  regardless  of  volatility 
considerations,  to  be  volume  distributed  above  100  microns.  The  model 
neglects  to  account  for  the  fact  that  the  specific  activity  begins  to 
fall  off  again  for  very  large  particle  sizes  (Tompkins,  1970:  382). 

The  resulting  activity  size  distribution  will  be  biased  with  the  very 
large  particles  carrying  too  much  activity  while  the  smaller  particles 
will  carry  too  little. 

Model  Description.  In  FrelHng  's  Radial  Power  Law  model,  the 
particle  size  distribution  was  assumed  to  be  log-normal.  Tompkins' 
modification  generalized  the  particle  size  distribution  to  allow  for 
arbitrary  distributions  to  be  entered  In  tabular  form.  Let 

dj^  «  geometric  mean  diameter  of  kth  particle-size 
,  class, 

b^  *  ^r.  refractory  fraction 

for  mass  chain  1, 

M|^  -  number  of  particles  In  class  k, 

Fy  >  total  number  of  fissions  In  all  size  classes, 
s  number  of  atoms  In  mass  chain  1  In  size  class  k, 
s  fission  yield  of  1th  mass  chain. 


iand 

r  ■  Index  for  a  refractory  chain. 

The  number  of  atoms  of  mass  chain  1  In  size  class  k  Is 


For  a  refractory  chain,  (b^  »  1),  Equation  (56)  becomes 

W-fT’r'H'V 

where  " 

v^k'  ■  Vk/£  "k^k 

/ 

Is  the  mass  fraction  In  the  kth  size  class.  The  fractionation  ratio 
Is  then 
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Equation  (61)  becomes 


I! 


It  then  follows  from  Equations  (56),  (60),  and  (64)  that 


"l<V  •  * 


Now  Tompkins' modification  assumes  a  two-component  syste^.  One  part 
obeys  the  Radial  Power  Law,  but  the  other  part  has  a  constant  specific 
activity.  Let 


»  fraction  of  atoms  In  Ivn  mass  chain  that  obeys 
the  Radial  Power  Law,  and 

«  fraction  of  atoms  In  1th  mass  chain  that  exhibits 
.  a  constant  specific  activity 


Equation  (65)  then  can  be  written 


;,((!„)  .  •  S,) 


The  crossover  point  between  the  Radial  Power  Law  and  the  constant 
distribution  yields 
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but  since 


or 


(68) 
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R^U  ♦  ^i^^cross^ 


(69) 


then 


1 


1-1 

1  *  ^1°cr^« 


(70) 


and 


bi-l 

^1*^cross 


1  ♦  E^O 


'1-1 


1  cross 


(71) 


Then  finally  we  have 


Nl(d,) 


®1-1 

1  *  E.D 

1  cross 


A  d| 


'1-1 


♦  D 


'1-1 


(72) 


11-35 


The  model  described  above  Is  termed  the  F-T  model  for  purposes  of 
comparison  with  the  G-X  model  to  be  developed  In  the  next  chapter. 

The  F-T  model  Is  considered  to  be  the  state  of  the  art  In  fallout 
modeling  by  present  researchers  In  this  area.  It  was  for  this  reason 
that  It  was  selected  for  comparison  with  the  G-X  model.  The  next 
section  will  describe  the  most  recent  effort  to  model  fractionation 
before  the  present  research. 

Diffusion  Model 

In  1966,  Korts  and  Norman  at  General  atomic  reconsidered  Miller's 
fractionation  model  to  Include  the  glassy  nature  of  fallout  particles 
from  a  general  land  surface  detonation.  The  occurrence  of  glass 
formation  suggested  the  use  of  temperature  variable  condensed  state 
diffusion  coefficients  In  place  of  the  temperature  switched 
equilibrium  model  proposed  by  Miller. 

When  describing  fission  product  sorption  according  to  a 
dlffuslvlty-condensatlon  model,  fallout  formation  Is  assumed  to  be 
governed  by  equilibria  established  at  the  gas-surface  Interfaces,  the 
rate  of  cooling  of  the  cloud,  and  the  rate  of  diffusion  of  adsorbed 
fission  products  Into  the  central  portions  of  the  fallout  particles 
(Norman,  1966;  12-54). 

The  assumption  that  soil  solidification  takes  place  at  a  fixed 
temperature,  as  Is  done  In  the  Miller  model.  Is  Inadequate  for 
silicate  soil  fallout  both  above  and  below  any  reasonable  "freezing 
temperature"  (Norman,  1966).  First,  "molten"  silicates  1n  the  region 
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Just  above  their  "freezing  temperatures*  are  generally  viscous  liquids 
which  have  low  mobilities  (diffusivltles)  of  their  component  species. 
Therefore,  above  the  "freezing  temperature,"  liquid  silicate  fallout 
particles  won't  be  uniformly  loaded  with  fission  products  but  will 
show  a  considerable  radial  concentration  gradient.  Additionally, 
molten  silicates  tend  to  form  glasses  on  cooling.  A  glass  Is  Just  a 
supercooled  liquid  whose  viscosity  has  become  so  high  that  the 
tendency  toward  crystallization  has  become  essentially  negligible. 
Norman  and  Vilnchell  suggest  that  diffusivlty  in  a  glass  can,  as  a 
first  order  approximation,  be  treated  as  an  extension  of  diffusivlty 
In  the  corresponding  molten  silicate  (Norman,  1966).  The  system 
behaves,  then,  as  If  no  phase  transition  has  occurred.  Thus  not  only 
Is  the  assumption  of  homogeneity  above  a  "freezing  temperature" 
inadequate,  but  the  assumption  of  strictly  surface  adsorption 
occurring  after  the  "freezing  temperature"  has  been  reached  Is  also 
questionable.  In  fact,  the  very  smallest  •allout  particles  will 
probably  absorb  fission  product  oxides  essentially  homogeneously  at 
temperatures  lower  than  the  "freezing  temperature"  (Norman,  1966). 
Using  temperature-dependent  diffusivltles  Instead  of  the  "freezing 
temperature"  model  Is  a  logical  extension  of  the  Miller  model. 

Sorption  of  fission  products  occurs  In  this  model  by  allowing  a 
homogeneous  gas  phase  to  equilibrate  with  the  surfaces  of  all  fallout 
particles.  The  adsorbed  fission  products  are  then  allowed  to  diffuse 
Into  the  fallout  particles.  All  of  this  occurs  In  the  time- 
temperature  field  associated  with  a  nuclear  detonation  (Norman,  1966). 


Concentration-Independent  diffusion  in  spherical  particles  as 
presented  by  Carslaw  and  Jaeger  (Carslaw,  1959)  is  assumed  to  govern 
the  diffusion  of  the  condensed  fission  products  into  the  carrier 
material.  The  diffusion  coefficient,  D,  is  related  to  the  average 
concentration,  and  the  surface  concentration,  Cg,  through  the 
following  equation: 


6  • 
I 

V  _  . 


1 

I? 


(73) 


/here  R  is  the  radius  of  the  sphe’^'cal  particle  and  t  is  the  diffusion 
time.  Both  the  Henry's  Law  constant  and  the  diffusion  coefficient  are 
assumed  to  follow  Arrhenius  and  Clausius-Clapeyron  temperature 
dependencies,  respectively  (Norman,  1966:  15). 

The  cooling  rate  as  a  function  of  yield  is  given  by  the  following 
cooling  rate  scaling  equation  given  by  Freiling  (Freiling,  1965): 

-  H  S8  3  X  10^^  W"°*^ 

(74) 

where  -dT/dt  is  the  cooling  rate  (degrees  Kelvin/sec)  and  H  is  the 
energy  yield  (Kt). 

The  mathematical  model  Involves  time-temperature  stepping, 

Henry's  Law  constants,  diffusion  coefficients,  a  detonation  model,  and 
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•  MSS  balance  equation  relating  the  total  amount  of  a  nuclide  In  the 
cloud  to  the  amounts  of  this  nuclide  In  the  gas  phase  and  In  the 
particles  of  various  sizes  through  an  analytically  determinable 
surface  concentration  of  the  nuclide.  By  stepping  the  temperature  and 
allowing  the  nuclides  to  decay*  a  similar  but  more  complicated 
detenel nation  of  an  Increment  to  the  surface  concentration  can  be 
Mde.  This  process  Is  then  repeated  until  the  diffusion  process 
essentially  ceases.  This  process  can  be  carried  out  through  each 
nuclide  chain  In  order  to  account  for  transmutation  effects.  In  this 
way*  it  is  possible  to  make  calculations  describing  the  chemlstrv  of 
fallout  formation  according  to  a  condensation-diffusion  mode  'Norman, 
1966:  21). 

Norman  achieved  some  promising  results  with  this  model  (Norman* 
1968).  According  to  Norman,  the  method  was  going  to  be  tested  for  use 
In  the  official  DoD  fallout  code*  OELFIC*  but  apparently  this  was 
never  done  (Norman,  1983).  It  was  because  of  the  very  encouraging 
results  that  Norman  reported  that  a  condensation-diffusion  approach 
was  chosen  as  a  basis, for  the  S-X  model  developed  In  the  next  chapter. 


III.  The  6-X  Diffusion  Model 


Introduction 

There  Is  a  considerable  amount  of  data  to  suggest  that  the 
following  description  of  the  fallout  formation  process  applies  to  a 
surface  burst:  Soil  particles  are  entrained  by  the  rising  fireball 
and  a  small  percentage  of  the  total  Is  completely  vaporized.  As 
condensation  begins,  very  small  particles  with  fission  products  fused 
within  them  are  created.  Many  of  these  coalesce  with  other  larger 
melted  and  unmelted  particles  along  with  the  vapors  of  the  more 
volatile  fission  products.  Those  soil  particles  which  are  not 
vaporized  generally  remain  In  the  fireball  for  a  shorter  period  of 
time.  In  addition,  they  constitute  the  larger  fallout  particles.  The 
unmelted  particles  enter  the  fireball  at  later  times  and  leave  the 
cloud  at  earlier  times.  Thus  they  receive  only  a  limited  amornt  of 
heating  and  consequently  the  fission  products  they  contain  are  fused 
into  the  outer  layer  or  attached  to  the  surfaces  only  {Miller,  1967: 

3,  Heft,  1970:  255-256,  Tompkins,  1970:  388,  Norment,  1966;  23-28). 

There  are  several  sources  of  Information  to  validate  this  picture 
of  the  early  fireball:  First  Is  the  separation  of  the  fallout  samples 
Into  two  distinct  types  of  particles:  glassy  and  crystalline  (hence 
the  name  G-X  model).  The  glassy  particles  are  those  which  were 
exposed  to  high  temperatures  and  either  fully  melted  or  vaporized. 

The  crystalline  particles  are  those  which  never  vaporized  or  melted  or 
at  most  only  slightly  melted.  Second,  the  two  types  of  particles 
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exhibit  very  distinct  behavior  with  respect  to  specific  activity.  The 
glassy  particles  which  were  exposed  to  higher  teriperatures  and  present 
In  Vie  cloud  earlier  both  absorb  refractory  fission  products  and 
adsorb  volatile  fission  products.  The  crystalline  particles  were  not 
present  In  the  active  region  of  the  cloud  at  early  time  and  so  absorb 
the  more  refractory  fission  products  to  a  much  smaller  extent.  This 
effect  Is  offset  some  by  the  tendency  of  the  smaller  particles  to 
agglomerate,  but  this  will  be  neglected.  This  model  of  behavior  Is 
supported  by  specific  activity  data.  Glassy  particles  In  general  show 
much  higher  specific  activity  than  crystalline  particles  of  the  same 
size  (Pascual,  1967:  8,  Mackin,  1958;  17,  Crocker,  1965:  5). 

In  addition  to  this  evidence  for  two  separate  types  of 
fractionation  behavior,  there  Is  evidence  that  the  simple  temperature 
switch  used  by  Frelling  and  Miller  Is  not  correct.  This  was  the 
premise  for  the  diffusion  studies  by  Norman  In  his  treatment  of 
fractionation  (Norman,  1966-1971).  Since  the  spherical  particles  are 
glasses  which  by  definition  are  supercooled  liquids,  they  exhibit  no 
well-defined  melting  point.  In  defense  of  his  diffusion  hypothesis, 
Norman  took  thin  sections  of  glassy  fallout  particles  and  obtained 
concentration  profiles  In  agreement  with  his  model  (Norman,  1967: 

219).  Figures  8,  9,  and  10  from  a  report  by  Miller  show  additional 
radiographic  evidence  for  diffusion  of  radioactivity  Into  fallout 
particles  (Miller,  1964:  22-31).  Additional  evidence  for  the 
diffusion  approach  from  specific  activity  measurements  on  ground 
bursts  Is  given  by  Nathans  (Nathans,  1970). 
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Figure  10.  Radiograph  of  an  Irregularly  Shaped 
Crystalline  Fallout  Particle  (Miller,  1963) 

in-5 


One  important  question  in  developing  the  diffusion  model  is 
whether  gas-phase  diffusion  of  fission  products  to  the  surfaces  of 
particles,  surface  attachment  coefficients,  or  diffusion  within  the 
particles  is  the  rate  controlling  mechanism.  In  a  working  group  at 
the  April  1966  Fallout  Phenomena  Symposium  (Proceedings,  1966:  456), 
Freiling  reported  that  work  done  by  U.S.  Naval  Radiological  Defense 
Laboratory  (NRDL)  on  calcium  ferrite,  •d0-Al202-Si02,  clay,  and 
sodium  oxide-silica  indicated  that  condensed  state  diffusion  was  the 
rate  determining  mechanism  for  all  but  the  calcium  ferrite  samples 
where  gas-phase  diffusion  was  found  to  be  rate  determining. 

Method  Overview 

For  the  G-X  model,  the  particles  are  assumed  to  be 
non-agglomerating,  non-convecting,  spherical  silicate  glassy  particles- 
of  various  sizes  which  are  distributed  uniformly  throughout  a  cloud 
which  is  at  a  uniform  temperature  and  has  a  uniform  concentration  of 
fission  products  and  air.  Calculations  proceed  on  one  decay  chain  at 
a  time_beginning  at  the  highest  temperature  at  which  diffusion  is 
considered  to  be  rate  determining.  At  this  point  the  isotopes  are 
assumed  to  surface  condense  according  to  Henry's  Law  and  then  to 
diffuse  into  the  fallout  particles.  Then,  the  isotopes  are  allowed  to 
decay,  and  the  time  and  temperature  are  incremented.  The  process  is 
repeated  until  the  temperature  is  close  to  the  solidification 
temperature  for  the  glass.  At  this  point,  the  distribution  of  glassy 
particles  is  augmented  with  the  distribution  of  crystalline 
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pTlfcUs.  The  proems  Is  conttmieP  ustil  the  diffusion  process  h,s 
esser.tUlly  temlneted.  The  techntques  appHed  fn  the  particle 
activity  module  here  «re  Influenced  by  the  «rk  of  Korts  and  Noman. 
The  propram  Is  on  file  at  the  Mr  Feme  Institute  of  Technolopy  and  Is 
fully  consented.  It  Is  deslpned  to  be  run  rith  2SO.OOO  dO-blt  «ords 
of  memory  In  three  overlays  on  a  COC  Cyber  750. 

Method 

hn  overview  of  the  method  has  been  presented.  This  section 
presents  a  derivation  of  the  Important  epuatlons  used  In  the  6-X  model 
For  an  assumed  Ideal  solution,  the  vapor  pressure  of  fission 
product  1  in  the  pas  around  a  llpuld  soil  particle  Is  piven  by 
Raoul t’s  Law  as 


where 

s  the  mole  fraction  of  1 
>  vapor  pressure  If  X.j  *  1,  that  Is  the 

vapor  pressure  which  would  be  found  above  pure 
fission  product  1  at  the  given  conditions. 

In  real  solutions,  there  are  often  departures  from  this  Ideal 
behavior.  Henry’s  Law  of  dilute  solutions  replaces  the  vapor 


pressure  of  the  substance  in  pure  form,  (the  proportionality 
constant),  with  a  different  constant.  This  is  called  the  Henry's  Law 
constant.  The  Henry's  Law  constant  is  the  slope  of  a  line  tangent  to 
the  observed  behavior  of  the  partial  pressure  of  fission  product  i  as 
a  function  of  at  very  low  concentrations  of  substance  i  (see 
Figure  4). 

^  *  *i  •‘i  (76) 

Substituting  for  P^  from  the  ideal  gas  law  we  have 

nJ^.X^k^  (771 

where 

R  «  gas  constant 
T  a  absolute  temperature 
\r a  volume 

or 

M^V 

•"I  ’  "I  ”i  *  ITT  ’^i  ''i  (78) 

where  , 
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•  MSS  of  fission  product  1  In  the  gas  phase 
■  Molecular  weight  of  1 


®son  * 

«  the  mass  of  1  In  solution 
a  the  molecular  weight  of  the  soil 

Because  of  the  very  low  concentration  of  fission  product  1  1n  the 
soil,  the  mole  fraction  Is  given  approximately  by 


(80) 


where  Is  the  concentration  of  1  In  the  particles. 

As  a  matter  of  convenience,  the  Henry’s  Law  constants  were 
measured  in  concentration  units  rather  than  In  molar  units.  That  is 


(81) 
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Thus  the  amount  of  fissfon  product  1  fn  the  gas  phase  at  equillbr-futn  Is 


- RT - 


(82) 


In  the  actual  situation,  diffusion  of  fission  product  atoms  to  the 
surfaces  of  the  particles  Is  generally  much  faster  than  diffusion  Into 
the  bulk  phase  of  the  particles.  In  this  case,  the  rates  of  uptake  of 
fission  products  are  controlled  by  the  diffusion  rates  In  the  bulk 
phase.  For  this  reason,  diffusion  through  the  gas  Is  taken  to  be 
Infinitely  fast.  In  addition,  because  the  cloud  Is  cooling  and 
because  the  fission  products  are  decaying,  equilibrium  cannot  be 
established  completely.  Instead,  the  gaseous  fission  products  are 
assumed  to  be  In  equilibrium  with  an  Infinitesimally  thin  surface 
layer  on  the  particles.  The  amount  of  fission  product  1  In  the  gas 
and  liquid  phases  Is  determined  by  a  nass  balance  as  follows:  At  the 
highest  temperature  at  which  flsslr-n  product  sorption  by  bulk 
diffusion  Is  determined  to  be  rate  detennining,  the  surface 
concentration  of  fission  product  1  Is  found  from 


C^l  k^j 

rt; 


m 


11 


where 


(83) 
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■  surface  layer  concentration  of  1  in  the  soil 
at  the  1st  time  and  temperature 
m  Henry's  Law  constant  1  at  1st  time  and 
temperature 

Vj  ■  volume  of  active  region  at  the  first  time 
Tj  «  temperature  at  the  1st  time 
m^^  s  mass  of  fp  which  has  diffused  into  the  soil 
particles 

Y^l  a  yield  of  fission  product  i  at  the  1st  time 

Th*  first  tern  on  the  right  of  Equation  (83)  is  the  amount  of  fission 
product  i  in  the  gas  phase.  The  second  term  on  the  right  must  be 
computed  by  solving  the  diffusion  equation  for  each  particle  size 
class  considered,  weighting  the  result  by  the  mass  of  soil  in  the  size 
class  and  summing  over  all  size  classes: 


P  (s) 

■"il  *  I  ^il  '^s 

s=l 


(84) 


where  is  the  average  concentration  in  the  particles  during  time 
period  one.  The  dif'^usion  equation  for  spheres  with  only  radial 
diffusion  is 

(85) 


On  letting 


u  •  C^r 


(86) 


Equation  (85)  becomes 


<i1th  boundary  conditions 

u  .  «  ^(t),  r  n  a,  t  >  0 
u  •  0.  t  -  0,  0  <  r  <  a 


Crank  gives  the  solution  In  this  case  as  (Crank,  1975) 

C  ,,.,).  f?li  (.1)*^  exp  (.D,kVt/a^^ 

-  ■  k.l 

'  (88) 

kTrsln(^) J  exp(0^kVx/a^)d(  x)dx 
x»0 

,Mch  1,  .  Volterr.  irtegr.!  .quatlon  of  the  seco-Kl  kind.  Sfoto  <(t) 
is  not  knovm  in  .<han«.  nn  U.r«1*e  technK^io  -ould  bo  roqoirod  to 
solve  this  equation.  For  this  reason,  the  solution  of  the  diffusion 
equation  Is  detemlned  by  a  superposition  technique.  For  the  first 


o 


;  ("> 


ttme  period,  the  surface  concentration  Is  held  constant,  and  the 
amount  of  fission  product  1  which  has  diffused  Into  a  particle  of  size 
s  during  the  time  period  Is  determined  by  solving 


n 


u«0,  r>0,  t>0 
u  -  a  r  .  a.  t  >  0 
uaO,  t>0,  0<r<a 

The  solution  to  this  equation  Is 


C(r.t),.  X  -)i 

I  k.l  *  ‘ 


or 


C(r,t)  » 


and  the  average  concentration  Is  given  by 


(89) 


(90) 


C^,(t)  »  ^  ^  ^  exp{-D^k  ir  t/a  )| 

'IT  (C 


ic«l 


(91) 
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Later  time  periods  are  treated  Independently,  and  the  results  are 
superimposed  to  give  the  final  result.  Returning  to  the  first  time 
period,  we  have  the  following  mass  balance: 


'11 


TJTT 


ViM^ 


♦  C 


11 


P 

T 

S..^ 


^11  “s 


(92) 


where  m^  Is  the  mass  of  soil  In  particle  size  class  s  and  p  Is  the 
number  of  particle  size  classes.  The  surface  concentration  Is  assumed 
to  be  Independent  of  particle  size  and  type.  At  this  point,  time  and 
temperature  are  Incremented.  To  find  the  surface  concentration  at  the 
second  temperature,  the  mass  balance  Is  given  by 


'12 


Cl  2  Vj  k^2 

- - 


£.(s) 


*  ^11  i^aJi;2  w^ 


*  C(1.1JJ  [l-exp(-x^  jAtj]  I  (93) 

Ssl 


♦  |C^2"*^il*’'P^"^f^^l^"^(1-l)l-^"*’'P^"^(1.1)*^l^^i  I  ^I^Lz  ''s 

Ssl 

The  first  terra  on  the  right  side  Is  the  amount  of  fission  product  1 
which  Is  still  In  the  gas  phase  at  tg.  The  second  term  Is  the 
amount  of  1  which  had  condensed  at  tj  after  being  allowed  to  decay 
for  4tj  and  to  diffuse  during  periods  atj  and  At2  (the  diffusion 
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coefficients  will  differ  for  the  two  periods).  The  third  term  is  the 
amount  of  the  immediate  precursor  of  Isotope  i  which  decayed  into  i 
and  was  diffusing  as  isotope  i-1  during  time  interval  at^  and  as 
Isotope  i  during  atg.  The  fourth  term  is  the  amount  of  isotope  1 
which  has  diffused  during  the  second  period  as  a  result  of  the  surface 
concentration  perturbation  caused  by  the  time-temperature>volume 
perturbation. 

Since  the  equations  at  succeeding  times  become  rather  involved, 
further  analysis  will  be  restricted  to  the  Instance  of  a  mass  chain 
containing  only  a  single  isotope.  Further,  the  discu^cion  will  be 
limited  to  the  case  of  a  single  particle  size  class  with  Wj^  «  1. 

These  restrictions  will  make  the  general  case  more  tractable.  In  this 
case.  Equation  (92)  becomes 


CjGj 


•  'l‘ll 


and  Equation  (93)  becomes 


Y 


C2G2  ♦  Cjaj2  *  ^''2‘’'1^^22 


(94) 


(95) 


or 


Y  »  C,6,  ♦  C 


2®2  *  V22  ^1^^12"*22^ 


(96) 


where 
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' 

/ 

o 

®M  ‘  ■'■  ■RTjT' 

At  the  third  time  period,  the  mass  balance  Is 

1 

Y  «  C3G3+CjAj3+(C2'Cj)a23  *  I!C3-(C2-C^)“Cj]a33 

(97) 

or 

-  •  - 

Y  .  C3G3+Cjai3*{C2-q)A23  ♦  (C3-C2)&33 

(98) 

At  the  fourth  time,  the  expression  becomes 

■<:) 

(99) 

♦  (€^-€3)^44 

At  the  Nth  step,  the  expression  Is 

t 

N 

k.2 

(100) 

or 

.. 

N-1 

Y  «  C^(V*Nn'  *  2,  ^k^*lcN'^lc+l,N^ 
kal 

(101) 

-  ■  : 
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In  order  to  show  the  effect  of  adding  crystal  particles  at  late 
time,  the  situation  must  be  made  slightly  more  complex.  Before  the 
crystalline  particles \;re  introduced,  the  expression  for  the  Nth  time 
step  is 


glass 


N-1 

iCsl 


z 

*glass 


kN 


)w 

Vl.N'''s 


(102) 


The  mass  balance  at  the  time  step  in  which  crystal  particles  are 
introduced  is 


Y  -  C  G  +0  V  A  *Y  C 
X  "x  X  ®xx  s  ^  ^ 


glass 


Z 

k>l  Is 


S/.(s)  (s)  V 


glass 


(103)- 


*'xZ 

*xtal 


^xx  ’^s 


‘If'-s-Z  Si 


all  s 


k>I  k 


glass 


(104) 


For  later  time  steps  (after  crystal  particle  entry),  the  expression 
becomes 


Ul-17 


''■‘'nS*  ItiZ  *s 

all  s 


M-1 


l*glass 


k«I  :s 


M-1 

*  Z'k 

ksx 


*xtal 


or  upon  rearrangement 


Y  •  Vm  *  ^mZ  ’^s 

all  s 


n-1 

•Z 

k»l 


‘^kM  ^kn,M''^s 

*gTass 


♦Z  <^42;  '4M’-4‘i,H»'s 

ksx  (all  s 


(105) 


(106) 


Each  term  (n  the  last  two  sums  In  Equation  (106)  contains  a  surface 

concentration  expression  times  the  difference  of  two  diffusion  terms. 

Each  of  these  must  be  modified  by  a  decay  term.  In  order  tc 

Illustrate  how  decay  Is  handled,  take  the  case  of  an  Isotope  with 

three  precursors  at  the  jth  time  step.  Figure  11  shows  the  ways  In 

which  growth  and  decay  can  occur  In  the  Instance  of  mass  chain  89.  In 
89 

the  case  of  Kr,  there  are  three  possible  paths  for  It  to 

89 

contribute  to  production  of  Rb  after  three  time  Increments.  The 
decay  coefficients  for  these  three  paths  are 
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Figure  11.  Decay  Paths  in  f'lass  Chain  89  for 
the  Diffusion  Model 
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(107) 


{ l-exp(  j )  1  exp{  -  a  13) ) 

{exp(-x^r^ti)||l-exp{-XnrAt2)|{exp(-Xj^jjAt3)|  ^^08) 

expC-Xj^^(Atj+At2)]|l-exp(-X|^^At3){  (109) 

Recall  that  each  decay  expression  is  associated  with  the  difference  of 
two  diffusion  terms  (see  Equation  (102;).  The  first  diffusion  temi 
represents  diffusion  from  tj^  to  the  current  step.  The  second  term 
represents  diffusion  from  tj^+j  to  the  current  time  step.  It  is 
seen,  then,  that  the  first  diffusion  term  corresponding  to  the  decay 
expression  in  equation  (107), is  for  diffusion  as  Kr  during  time  step  1 
and  as  Rb  from  t2  to  t^.  The  second  diffusion  term  is  for 
diffusion  as  from  t2  to  t^.  Similarly,  the  first  diffusion 
term  associated  with  Equation  (108)  is  for  diffusion  as  Kr  from  tj 
to  tg  and  as  Rb  from  t3  to  t^.  The  second  term  is  then  for 
diffusion  as  Kr  during  At2  and  as  Rb  from  t3  to  t^.  And, 
finally,  the  first  diffusion  term  associated  with  Equation  (109)  is 
for  diffusion  as  Kr  from  t^  to  t3  and  as  Rb  during  At^,  while  the 
second  is  for  diffusion  as  Kr  from  t2  to  t3  and  as  Rb  during 
At^.  Each  of  these  diffusion  terms  is  just  the  solution  of  the 
diffusion  equation  using'  the  appropriate  argument.  For  example,  the 
first  d’lffusion  term  associated  with  Equation  (107)  was  for  diffusion 


n:-2o 


V  , 


as  Kr  during  time  step  1  and  as  Rb  from  t2  to  t^.  The  appropriate 
argument  for  the  diffusion  equation  Would  be 

.  •  0lCr.l*h^°Rb.2'‘V°Rb.3''V'^Rb.4^^4  (HO) 


Equation  (91)  gives  the  relationships  between  the  average 
concentration  and  the  surface  concentration.  The  derivation  of  the 
general  case  for  the  Hth  time  step  is  analogous  to  the  derivative  for 
Equation  (102).  For  times  prior  to  the  Introduction  of  crystal  soil 
particles,  the  average  concentration  is 


S  .  ksl 


glass 


'glass 


(111) 


Expressions  analogous  to  Equations  (103)  and  U04)  can  also  be 
similarly  derived. 

Particle  Size  Distribution 

The  lognormal  distribution  given  by 


/i  ,  l,lnx-ax 


(112) 


where 


a  cunnjlatlve  mass  fract<or 

fli 

K  >  a  normalization  constant 
a  •  In  (median  diameter) 

B  a  logarithmic  standard  deviation  of  the  distribution 
d  a  particle  diameter 

has  been  traditionally  employed  In  fallout  calculations.  The 
experimental  basis  for  such  a  choice  Is  found  In  the  analysis  of  cloud 
and  fallout  samples  obtained  several  hours  after  burst  (Heft,  1970, 
Norment,  1966:  28-29).  Based  on  an  analysis  of  data  from  the  Small 
Boy  event  In  1962,  a  lognormal  distribution  of  particle  size  with 
median  diameter  of  .407  microns  and  a  standard  deviation  of  In  (4.0) 
was  selected  as  the  default  distribution  for  the  DELFIC  code 
(McDonald,  1974:  60,  Kalonel,  1974).  In  this  distribution,  50  percent 
of  the  mass  Is  associated  with  particles  smaller  than  130  microns  In 
diameter.  Later  studies  of  the  soil  and  debris  ejected  from  the 
craters  of  various  tests  where  particle  sizes  ranged  from  millimeters 
to  meters  Indicated  a  power  law  behavior  (Layson,  undated): 

M(d)  .  k'  d'P  (113) 

where  p  Is  approximately  0.5  for  various  soils  and  rock.  In  1966, 
russell  suggested  a  truncated  power  law  based  on  evidence  that  a  true 
power  law  was  unrealistic  (Bridgman,  1983  d).  A  year  later,  Frell.ing 
published  a  study  which  concluded  that  the  differences  between  the 
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lognormal  distribution  and  the  truncated  power  law  were 
Insignificant.  He  further  concluded  "that  the  lognonnal  distribution 
had  the  esthetic  advantage  of  an  observational 1y  confirmed  theoretical 
basis*  (Frelling,  1966).  It  seems,  however,  that  requiring  the  active 
particles  which  obtained  their  activity  from  the  fireball  (as  opposed 
to  activation  products)  to  fit  the  distribution  of  very  large 
particles  ejected  from  the  crater  Is  erroneous.  At  any  rate,  one  can 
attempt  to  bridge  the  gap  between  these  two  positions  In  spite  of 
Frelllng's  findings  by  using  a  hybrid  distribution  as  suggested  by 
McGahan  (McGahan,  1974:  61).  McGahan  has  done  calculations  with 
lognormal,  pure  (not  truncated)  power  law,  and  a  hybrid  distribution 
and  discovered  some  differences  In  predicted  dose  rates.  But  If  one 
substitutes  the  more  reasonable  truncated  power  law  distribution, 
there  Is  little  difference.  As  a  result,  for  the  purpose  of  this 
study,  the  preferred  particle  size  distribution  Is  the  lognonnal. 

This  Is  especially  so  In  view  of  the  fact  that  this  research  Is 
compared  primarily  with  shots  Small  Boy  and  Johnny  Boy  %hich  were  of 
similar  yield  and  shot  conditions  and  since  a  lognonnal  was  deemed  to 
best  fit  the  Small  Boy  data. 

Partition  of  Soil 

The  next  Important  question  to  consider  Is  that  of  the  partition 
of  the  soil  burden  of  the  cloud  Into  the  various  phases  at  early 
time.  This  Is  a  question  which  would  best  be  answered  by  a  combined 
theoretical  and  computational  effort  that  would  start  from  first 


principles  and  follow  the  numerous  complex  processes  from  shot  time  to 
early  cloud  development.  Because  of  the  enormous  complexities 
Involved,  however,  no  such  calculation  has  been  attempted,  and  It 
remains  an  Important  area  for  further  research.  On  the  other  hand, 
bounds  can  be  placed  on  the  problem.  The  total  soil  burden  of  the 
cloud  can  be  estimated  from  a  number  of  sources:  (1)  energy 
partition,  (2)  hydrodynamics  calculations,  (3)  cloud  bouyancy,  (4) 
back  extrapolation,  (5)  specific  activity,  (6)  crater  volume,  and  (7) 
fireball  volume.  Norment  and  others  studied  and  evaluated  all  of 
these  and  chose  methods  4-7.  The  scaling  functions  they  developed 
were  calibrated  with  data  from  the  Teapot  Ess  shot  which  was  a  1  kt 
underground  burst  '.-67  feet)  at  the  Nevada  Test  Site  In  1955 
(Hawthorne,  1979:  201). 

The  specific  partition  of  the  soil  Into  vaporized,  melted,  and 
unmelted  material  Is  more  difficult.  First  of  all,  the  cloud  Is 
anything  but  uniform  at  early  time.  It  Is  extremely  hot  In  the  region 
cf  the  toroidal  ring  about  which  dust  and  debris  are  being 
circulated.  The  temperature  falls  off  rapidly  with  Increasing 
distance  from  the  center  of  this  region.  It  Is  reasonable  to  expect 
that  the  fission  products  and  device  debris  would  be  found  In  the 
hottest  regions  of  the  toroid.  Similarly,  one  would  expect  to  find 
any  soil  which  had  been  Initially  vaporized  located  here.  As  the 
fireball  rises.  It  expands  and  thus  heats  and  contaminates  additional 
soil.  After  an  initial  pseudohydrostatic  cloud  rise  phase  (the  second 
temperature  maximum  occurs  In  this  phase)  the  fireball  begins  to  rise 


rapidly.  A  strong  updraft  Is  produced  In  the  wake  of  the  cloud,  and 
soil  dust  suspended  1n  the  air  Is  sucked  upward  In  the  stem  toward  the 
cloud  cap  (Figure  1). 

The  next  thing  to  consider  Is  the  soil  Itself.  For  continental 
United  States  soils,  the  characteristic  soil  types  are  the  common 
sandy  and  clay  soils.  In  the  case  of  sand,  the  chemical  composition 
Is  nearly  100  percent  silicon  dioxide  (SlOg).  For  the  more  abunda^'t 
clay  soils,  the  major  component  Is  still  silicon  dioxide,  but  there  Is 
a  mixture  of  AI2O,,  CaO,  Na2C,  Pe2®3' 

SbgO^  as  v^ell.  In  studies  of  thermodynamic  properties  and 
diffusivltle?  of  soils,  Homan  and  others  concentrated  on  a  mixture  of 
62  percent  SIO2,  15  percent  AKO^  and  23  percent  CaO,  although 
other  mixes  were  studied  too  (Homan,  1970,  Homan,  1966,  Homan, 

19C7,  Vlinchell,  19C7).  Since  It  was  detemlned  that  this  mix  was 
appropriate  for  Small  Boy  calculations  (Momen,  1968:  1),  this  Is  the 
soil  mix  which  was  selected  for  this  study.  The  major  constituents  of 
this  soil  are  silicon  dioxide  and  calcium  oxide.  The 
Cao-Al^Oj-SiOg  soil  "solidifies"  at  about  1620  K.  To  estimate 
the  amount  of  soil  material  which  Is  In  the  vapor  phase,  the  following 
observations  are  made:  (1)  High  temperatures  are  localized  to  the 
inside  of  the  vortex  toroid.  (2), This  region  contains  most  of  the 
fission  products.  (3)  In  surface  bursts,  much  of  the  fallout  Is 
composed  of  unvaporized  soil  onto  which  active  material  condenses. 

Fror.i  these  it  follows  that  the  amount  of  soil  vapor  In  the  cloud  will 
be  less  than  the  mass  of  active  fallout.  In  shots  Jdhnny  Boy  and 
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Small  Boy,  about  38  percent  of  the  gross  fallout  mass  Is  associated 
with  active  material.  Therefore,  It  follows  that  at  least  for  these 
shots,  less  than  38  percent  of  the  total  soil  burden  Is  vaporized 
(Norment,  1966).  Norment  et  a1.  made  estimates  of  the  energy 
available  for  heating  soil  by  differencing  the  thermal  energies 
radiated  by  air  and  surface  bursts  and  found  that  1.5  percent  of  the 
soil  burden  Is  vaporized  per  100  degree  excess  of  the  temperature  over 
the  soil  boiling  point  at  a  time  shortly  after  the  second  temperature 
maximum  (Norment,  1966:  28).  Tompkins,  et  al.  have  Inferred  from  the 
magnitude  of  specific  activity  of  fallout  from  ground  bursts  that  only 
about  a  tenth  of  the  soil  Is  vaporized  (Tompkins,  1968). 

Crystal -61  ass  Distributions 

In  examining  samples  of  Small  Boy  fallout  which  w?r*>  not  sieved, 
Pascual  noted  bimodallty  In  the  sample  collected  at  9200  feet  from 
ground  zero  (Pascual,  1967).  Besides  this  and  other  qualitative 
remarks  made  above  by  various  Investigators,  there  are  a  number  of 
reports- which  provide  some  quantitative  Information  on  the  partition 
of  the  particle  size  distribution  between  glassy  particles  and 
crystalline  particles.  Two  reports  which  looked  specifically  Into 
this  are  Pascual,  1967,  and  Mackin,  1958.  Pascual  provides  the  most 
relevant  data  since  It  uses  Small  Boy  data.  Unfortunately,  only  the 
larger  particle  size  classes  were  Investigated  In  the  report. 
Information  on  the  smaller  cl isses  had  to  be  inferred  from  knowledge 
of  the  overall  particle  size  d'- jtrlbutlon  and  other  Information. 


Table  1  gives  the  raw  data  used  to  determine  the  fraction  of  each  size 
class  which  Is  composed  of  glassy  particles  as  opposed  to  crystalline 
particles.  The  two  distributions  were  unfolded  by  first  assuming  an 
overall  lognormal  particle  size  distribution  of  median  .407  microns 
diameter  and  standard  deviation  of  In  (4.0).  Then  using  the 
Information  In  Table  1,  qualitative  remarks  such  as  that  spheres 
dominate  In  the  smaller  particle  size  classes  (Nathans,  1970:  367), 
and  specific  activity  data,  the  two  distributions  were  determined 
using  a  non-linear  search  techni.^t.'e  for  the  parameters  of  two 
lognormal  distributions  which  when  s'.nerimposed  would  yield  the 
closest  fit  to  the  observed  combined  d  tributlon.  The  resulting 
distributions  are  shown  In  Figure  12. 

Bias  In  the  Data 

There  was  some  concern  that  measurements  on  the  particle  size 
distribution  for  shot  Small  Boy  were  in  error  oue  to  the  sieving 
technique  used.  The  concern  was  that  in  the  process  of  separating  the 
fractions  using  sieves  of  different  mesh  sizes  that  small  particles 
adhering  to  the  surfaces  of  the  larger  particles  would  be  shaken  loose 
and  thus  bias  the  analysis.  Also  the  larger  particles  in  some  cases 
were  fragile  and  could  be  broken  during  size  separation.  Pascual 
looked  Into  this  source  of  potential  error  and  concluded  that  while 
there  was  measurable  loss  of  activity  from  the  larger  size  fractions, 
the  gamma  ray  spectra  were  not  noticeably  altered  (Pascual,  1967). 


ni-27 


Table  1 


Raw  Data  for  Mass  Size  Distribution  Decomposition 

Parti cl  * 

Size 

Particle 

Type 

Percent 

In  Class 

Specific 

Activity 

CPM/mg 

' 

Spheres 

35.7 

3617 

1400m 

Crystal s 

64.2 

1657 

■ 

• 

Spheres 

38.4 

5677 

•700-1400y 

Crystals 

61.6 

2780 

Spheres 

46.1 

6700 

350-700m 

Crystal  s 

_ _ 

53.9 

3240 

(Pascual,  1S67) 
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Miller  and  Yu  accounted  for  this  effect  when  determining  the  total 
specific  activity  as  a  function  of  particle  size  (Miller,  1967). 

Determination  of  Diffusion  Coefficients 

Norman  and  others  spent  considerable  effort  In  determining 
diffusion  constants  for  various  soil  types  (Norman,  1966-1971).  The 
soil  type  of  Interest  here  Is  the  CaO-Al^O^-SIOg  soil  which  Is 
chosen  as  representative  of  common  soils. 

The  Arrhenius  temperature  dependence  for  diffusion  coefficient 
may  be  written  as 


D 


®o  « 


.E*/RT 


(114) 


where  D^  and  E*  (activation  energy  kcal/mole)  are  found 
empirically.  Data  collected  by  Winchell  and  Norman  Indicated  a  close 
correlation  between  E*  and  D^  so  that  If  one  Is  able  to  make  an 
Independent  determination  of  one,  one  can  compute  the  other  with 


logioDo  =  -5.79  *  1.23  x  10'^  E* 


(115) 


This  relationship  Is  known  as  the  compensation  law  (Winchell,  1967: 
487).  As  a  means  of  estimating  D^,  one  very  successful  method  Is  to 
use  the  Ionic  radius,  r,  of  the  diffusing  Ionic  species.  This  applies 
strictly  only  In  the  case  of  monovalent  Ionic  species,  but  Winchell 
and  Norman  suggest  using 


p  ■  1«59  r 


(116) 

where  r  Is  the  Ionic  radius  for  the  most  probable  oxidation  state  In 
the  silicate  1*  high  temperature.  Table  2  lists  the  constants  as 
determined  using  these  procedures  for  all  of  the  significant  diffusing 
species.  The  data  on  Ionic  radii  are  taken  from  Weast  (Weast,  1974). 

Determination  of  Henry *s  Law  Constants 

Norman  et  al.  have  made  numerous  studies  In  order  to  determine 
Henry's  Law  constants  for  those  atomic  species  of  Interest  In  fallout 
research  (Norman,  1966-1971).  The  values  used  In  this  research  are 
listed  In  Table  3.  In  several  cases,  no  data  was  available.  These 
few  cases  were  for  less  Important  species  but  for  completeness, 
estimates  of  the  Arrhenius  coefficients  were  made  for  these  by  ' 
comparlrig  them  with  those  species  which  demonstrated  similar 
volatility.  These  cases  are  Indicated  In  the  table. 

Theoretical  Considerations 

While  many  phenomena  remain  to  be  Incorporated  Into  any 
fractionation  model,  one  which  should  be  given  attention  Is  the 
apparent  decrease  In  mixing  efficiency  between  the  hot  gases  In  the 
fireball  and  the  ioll  as  yield  Increases.  A  hydrodynamic  calculation 
would  be  required  to  determine  the  extent  of  this  effect.  Fractional 
condensation  does  not  appear  to  be  nearly  as  Important  In  ground 
surface  bursts  as  It  Is  In  air  or  tower  bursts. 
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Atom-fc  No 


OCl 


0C2 


27 

2.54 

14800 

28 

-4.43 

2400 

29 

2.19 

14200 

30 

-3.3 

4400 

31 

.65 

11400 

32 

-1.44 

7700 

33 

-.0237 

9900 

34 

-2.13 

6500 

35 

8.64 

36 

2.01 

!fi  11 

37 

4.27 

18200 

38 

2.01 

13900 

39 

-.618 

9200 

40 

3.16 

15900 

41 

3.77 

17000 

42 

.94 

12500 

43 

.46 

11100 

44 

-4.93 

1500 

45 

-4.68 

1900 

46 

3.23 

16000 

47 

6.56 

21900 

48 

4.83 

18900 

49 

3.57 

18300 

50 

2.49 

15700 

51 

2.1 

14100 

£2 

3.3 

15600 

53 

7.29 

22500 

54 

3.51 

16500 

55  . 

5.57 

20500 

56 

6.91 

22600 

57 

5.2 

19500 

58 

4.42 

18100 

59 

4.25 

17800 

60 

.712 

11500 

61 

.522 

11200 

62 

.338 

10900 

63 

4.68 

18600 

64 

.005 

10300 

65 

4.47 

18200 

66 

-.403 

9600 

Atomic  No. 

Boiling  Point 

HCl 

HC2 

27 

3173 

9.2 

33500 

28 

3173 

(9.2) 

(33500) 

29 

2907 

(9.7) 

(33500) 

30 

3000 

11.8 

23000 

31 

2976 

9.7 

35200 

32 

1764 

12.5 

23600 

33 

1010 

17800 

34 

1026 

7.0 

!  4200 

35 

331.8 

4.0 

i  5200 

36 

120.1 

4.0 

500 

37 

1656 

7.0 

1  18100 

38 

3497 

4.4 

26700 

39 

4695 

9.4 

,  5PFC0 

40 

4808 

6.7 

38200 

41 

3300 

8.5 

44800 

42 

1351 

7.6 

20400 

43 

583 

4.1 

7100 

44 

4505 

1.2 

2600 

45 

4149 

7.2 

,  18400 

46 

3436 

8.7 

24700 

47 

2451 

6.2 

14300 

48 

1831 

9.0 

19200 

49 

2123 

9.1 

30600 

50 

2247 

13.1 

32100 

51 

1832 

11.1 

29500 

52 

1534 

9.2 

11800 

53 

457.4 

3.8 

1000 

54 

165.9 

4.1 

700 

55 

1555 

6.6 

17300 

56 

•  3003 

4.8 

22500 

57 

4608 

7.3 

41800 

58 

4367 

7.0 

35400 

59 

4252 

5.7 

33800 

60 

4464 

7.6 

35700 

61 

4348 

7.3 

37400 

62 

4300 

6.7 

41400 

63 

4300 

(6.7) 

(41400) 

64 

4300 

(6.7) 

(41400) 

65 

4300 

(6.7) 

(41400) 

66 

4300 

(6.7) 

(41400) 

See  text  for  references.  Numbers  in  parentheses  are  estimated 


IV.  Results  and  D1scus!;Ion 


Calculations  were  made  with  the  standard  DELFIC  computer  program 
and  with  the  G-X  Model.  In  addition  to  these,  computations  to 
generate  more  detailed  data  on  mass  chains  89  and  95  were  made  with  a 
modification  of  the  G-X  code.  The  output  data  at  the  time  of  Interest 
Included  number  densities  for  the  nuclides  of  Interest  as  a  function 
of  particle  size,  specific  activities  for  those  nuclides  as  a  function 
of  particle  size,  total  yield  of  those  nuclides,  the  activity  size 
distribution,  and  the  dose  size  distribution.  These  raw  data  were 
than  taken  and  values  as  a  function  of  particle  size  were 
determined  for  each  method.  Log-log  plots  of  R^  gg  versus  Rgg  gg 
values  were  then  prepared  and  fit  with  linear  least  squares  to 
detemlne  the  correlation  slopes  (Bevington,  1969:99-102). 

In  the  plots  of  Rgg  gg  as  a  function  of  particle  size  (Figure 
13),  It  Is  Immediately  obvious  that  the  G-X  method  Is  superior  to  the 
pure  Frelling  method  up  to  about  100  microns,  but  below  about  30 
microns,  the  Frelling-Jompkins  (standard  DELFIC)  Model  more  closely 
matches  Small  Boy  data.  Also  at  la^e  particle  sizes,  the  G-X  method 
is  far  superior  to  both  methods.  In  terms  of  the  logarithmic 
correlations  which  are  popular  In  fractionation  analysis.  Tables  4  and 
5  list  the  correlation  slopes  for  each  of  the  methods  ard  compares 
them  to  Small  Boy  and  Johnny  Boy  data  respectively  (Frelling,  1968, 
crocker,  1965:  78)).  In  most  cases  the  agreement  between  the  G-X  . 
Method  and  the  data  Is  as  good  as  the  other  methods.  In  the  case  of 
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Figure  13.  Rgg  gg  versus  Particle  Size  for  Each  Method 
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Mass  Chain 
89 


Small  Boy 
1 


6- 


Freiling-Tomp. 

1 


90 

.73 

91 

.51 

95 

.0 

99 

.04 

.31 

.84 

132 

.90 

136 

.65 

137 

1.19 

140 

.5.': 

1 


.83 

- 

.64 

.40 

.0 

.0 

o 

« 

1 

.0 

.58 

.17 

.61 

.31 

1.39 

1.49 

.91 

1.45 

.34 

.73 

144 


03 


-.14 


01 


Table  5.  Logarithmic  Correlation  Slopes  -  Johnny  Boy 


Mass  Chain 

Small  Boy 

6-X 

Frelling-Tomp. 

89 

1 

1 

1 

90 

- 

- 

91 

.54 

.576 

.422 

95 

.0 

.0 

99 

.017 

-.028 

.006  . 

131 

1.10 

.459 

.169 

132 

1.08 

.491 

.313 

136 

.83 

.491 

1.48 

137 

1.11 

.916 

1.45 

140 

.61 

.304 

.736 

144 

1 

• 

o 

o 

-.052 

.022 

chains  131  and  13?  there  Is  considerable  Improvement  over  the 
Frelllng-Tompklns  model.  The  1og>1og  plots  for  both  Small  Boy  and 
Johnny  Soy  are  Included  in  Appendices  B  and  C,  respectively.  Figure 
14  shows  *he  total  specific  activity  for  the  various  methods  along 
with  th(?  measured  specific  activity  data  for  shot  Small  Boy  (Miller, 
1967:33).  Only  the  S-X  Model  follows  the  drop  In  specific  activity 
for  the  larger  particle  sizes.  All  methods  predict  an  Increase  In 
specific  activity  with  decreasing  particle  size  for  the  smaller 
particles.  This  fact  has  been  observed  by  many  Investigators,  but 
accurate  specific  activity  for  the  smaller  particle  sizes  Is  not 
available,  only  the  tendency  In  the  data.  Figure  15  shows  the  dose 
size  distributions  for  the  different  methods.  The  comparisons 
Indicate  that  the  Frelllng-Tompklns  model  underpredicts  the  down  field 
dose  rate  while  overpredicting  the  close  In  dose  rate.  As  mentioned 
In  the  introductory  chapter,  this  fact  will  have  Important 
consequences  for  lethality  calculations. 

Plots  of  specific  activity  for  specified  chains  are  also 
Included,  In  Apjpendix  ,0.  They  show  clearly  the  effect  of  treating  the 
crystalline  particles  separately  from  the  glassy  particles.  Perhaps 

89  Qc 

more  striking  are  the  plots  of  the  concentration  of  Sr  and  Zr 

for  each  type  of  particle  as  a  function  of  particle  size.  Figures  16 
89 

and  17,  Since  Sr  Is  a  volatile  species.  Its  distribution  on  the 
crystalline  particles  Is  very  similar  to  the  distribution  In  the 
glassy  particles.  For-  "'Zr,  on  the  other  hand,  there  Is  a  marked 
drop  In  the  concentration  In  the  crystalline  particles  for  the  larger 
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Figure  16.  Strontium  90  Concentration  as  a  Function  o' 
Particle  Size  for  Crystalline  Particles,  Glassy  Particles, 
and  the  Weighted  Average 
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particle  sizes.  The  effects  of  diffusion  are  evident  at  both  ends  of 

the  size  spectrum.  For  small  particles,  even  crystalline  particles 

permit  volume  loading.  For  large  particles,  diffusion  into  the 

particle  interiors  is  time  consuming  and  these  particles  cannot  load 
95 

uniformly  vdth  2r  even  for  the  glassy  particles. 

Figure  18  shows  the  total  spacific  activity  for  each  of  the 
particle  types.  The  crystalline  particles  show  a  clear  surface 
distribution  (decrease  in  specific  activity  with  increasing  particle 
size),  while  the  glassy  particles  are  surface  distributed  for  small 
particles  and,  very  large  particles  and  are  volume  distributed 
(constant  specific  activity)  in  the  intermedicite  size  range.  The  drop 
in  specific  activity  for  the  large  glass  particles  is  probably  due  to 
the  slov/  uptake  of  refractory  nuclides  by  the  spheres  due  to  large 
particle  size.  .  .  - 

Figures  19  and  20  show  a  comparison  of  the  observed  slopes  to  the 
calculated  slopes  for  Small  Boy  and  Johnny  Boy  respectively.  The 
error  magnitudes  are  plotted  as  histogram  data  in  figures  21  and  22. 
For  shot  Small  Boy  the  average  error  for  the  6-X  method  was  .245 
versus  .265  for  the  F-T  method.  For  shot  Johnny  Boy  the  average  error 
was  .276  for  the  G-X  method  versus  .341  for  the  F-T  method.  Thus  in 
both  shots  the  G-X  method  gave  better  resul cs. 
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V.  Conclusions  and  Recommendations 


Conclusions 

While  the  application  of  these  models  Is  restricted  to  ground 
surface  bursts  over  silicate  soils,  they  still  have  a  wide  range  of 
applications.  The  assumptions  male  In  developing  the  models  are 
considered  to  be  satisfactory  In  view  of  the  many  sources  of  error. 

As  It  Is,  the  assumptions  made  here  are  somewhat  less  restrictive  than 
those  for  the  previous  models,  although  there  Is  still  room  for 
Improvement.  The  success  of  the  £-X  Model  for  the  two  shots 
considered  Is  very  encouraging,  however.  While  the  S-V  model  (see 
Appendix  A)  shows  promise  for  faster  operational  type  codes,  the  G-X 
Model,  because  o*  the  overall  better  performance.  Is  the  preferred  one 
for  work  requiring  more  rigor,  especially  If  concentrations  of 
specific  Isotopes  are  required.  The  6-X  Model  certainly  has  the 
advantage  of  being  based  on  physics  rather  than  being  an  empirical  fit 
as  Is  the  radial  power  model. 

Recommendations 

There  are  a  number  of  areas  which  need  further  work  In  the 
fallout  area.  If  one  were  Intent  on  doing  the  problem  correctly,  all 
of  the  following  things  would  need  to  be  done: 

(1)  The  particle  size  distribution  needs  to  be  determined 
dynamically  from  early  time  In  the  fireball  to  Include  the  partition 
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between  solid,  llqi'ild,  end  gas.  E^>has1s  should  be  plateo  on 
tetenalnlng  the  extent  and  effect  of  agglomeration. 

(2)  The  degree  of  mixing  of  the  f1s«1on  products  with  the  soil 
carrier  as  a  function  of  time  and  Iccatlon  In  the  hot  toroid  should  be 
determined. 

(3)  Henry  ’s  Law  constants  a.id  diffusion  coefficients  need  to  be 
measured  for  a  variety  of  soil  types. 

m 

(4)  Interdiffusion  coefficients  for  the  diffusion  of  the  fission 
product  gases  through  the  air  to  the  surfaces  of  the  soil  particles 
should  be  determined. 

(51  Surface  attachment  coefficients  for  each  Important  compound 
need  to  be  measured. 

(6)  An  Investigation  of  fall  rates  for  different  particle  types 
needs  to  be  done. 

(7)  Work  on  the  effect  of  fractional  condensation  tnd  compound 
formation  by  fissicn  product  oxides  with  the  carrier  material  should 
be  Investigated. 

(•8)  Wore  Information  Is  needed  on  the  thermal  stability  of  the 
fission  product  oxides. 

(9)  And  finally,  work  needs  to  be  done  on  the  Influence  of  water 
on  the  reaction  balances. 

If,  on  the  other  hand,  one  Is  more  concerned  with  Improving  the 
faster  running,  so-called  “smear*  codes,  this  can  be  accomplished  by 
several  techniques:  One  Is  to  simply  take  the  S-Y  Model  of 
fractionation  and  determine  decay  rates  for  each  of  the  components 
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(surface  activity  and  volume  activity).  Alternatively,  one  cwuld 
follcw  only  those  chains  which  rake  the  largest  contribution  to  the 
surface  and  volume  components  of  the  activity  and  distribute  that 
activity  using  the  S-V  Model.  Then  the  result  could  be  normalized  so 
that  the  combined  total  activity  agrees  with  that  which  would  be 
computed  by  the  Uay-WIgner  formula.  Either  of  these  two  techniques 
would  allow  for  an  explicit  calculation  of  fractionation  effects  In  a 
code  such  as  the  Air  Force  Institute  of  Technology  model  (Bridgman, 
198:). 

Finally,  the  6>X  Model  could  be  used  In  a  parametric  study  of 
yield  and  fuel  type  (and  perhaps  other  parameters)  to  determine  a  new 
set  of  activity  size  distributions  for  SEER,  LASEER,  PROFET,  and  KDFOC 
which  are  all  derived  from  DELFIC  calculations.  These  codes  do  not  do 
explicit  particle  activity  calculations  but  rather  rely  on  fits  of 
OELFIC  calculations. 
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Appendix  A 


The  S-V  Model 

Introduction 

Because  of  the  ease  and  speed  with  which  a  model  developed  by  the 
Air  Force  Institute  of  Technology  (Bridgman,  1982)  can  be  used,  and 
because  the  Air  Force  Institute  of  Technology  (AFIT)  model  has  the 
potential  for  accounting  for  fractionation  In  an  explicit  manner,  the 
first  mode!  developed  as  part  of  this  research  was  one  which  would 
allow  direct  adaptation  for  use  with  the  AFIT  model.  The  new  model 
described  here  wHT  be  termed  the  S-V  model  because  It  reparates  the 
activity  size  distribution  Into  surface  and  volume  distributed 
components. 

Model  Description 

The  S-V  model  Is  a  modification  of  Miller's  method.  The  method 
Is  Implemented  by  first  computing  the  time  at  which  the  cloud  would 
reach  the  soil  solidification  temperature.  The  fission  product 
1nventOT7  Is  then  calculated  at  the  time  of  soil  solidification  using 
rigorous  computations  with  the  Bateman  equations.  After  that,  the 
amount  of  soil  In  the  active  region  of  the  cloud  Is  computed  along 
with  cloud  volume.  Next,  the  amount  of  material  which  will  dissolve 
Into  a  bulk  liquid  soil  at  that  temperature  Is  computed  using  Henr7's 
Law  (these  calculations  are  discussed  In  the  next  sub-section).  This 
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portion  of  the  fission  product  Inventory  Is  then  saved  temporarily  In 
a  large  array.  Those  fission  products  which  remained  In  the  gas  phase 
are  also  saved  In  another  array.  The  Bateman  equations  are  solved 
again  to  compute  the  decay  of  the  solid  phase  fission  products  In  the 
first  array  to  the  time  of  Interest,  taken  here  to  be  H  +  1  hour. 

Those  decay  products,  regardless  of  their  volatility,  are  assumed  to 
be  locked  Into  a  solid  matrix  and  are  distributed  uniformly  throughout 
the  volume  of  the  fallout  particles. 

Next  the  second  array,  representing  the  gas  phase  at  soil 
solidification.  Is  allowed  to  decay  to  H  *  1.  These  fission  products 
are  then  distributed  according  to  the  distribution  of  surface  area  for 
the  fallout  particles.  The  assumption  here  Is  that  the  particles  are 
now  solid  and  will  not  absorb  any  more  fission  products  but  may  adsorb 
them.  A  further  assumption  Is  that  no  particles  leave  the  cloud  until 
all  fission  products  have  been  absorbed.  A  refinement  of  the  method 
allows  the  last  assumption  to  be  relaxed  and  Is  discussed  In  the 
recommendations. 

D1 scusslon 

The  model,  while  It  Is  very  simplistic,  does  account  In  a  direct 
way  for  the  major  processes  Involved  In  particle  formation.  Its  major 
drawbacks  are  that  It  does  not  account  for  the  presence  of  crystalline 
particles  nor  for  the  fact  that  larger  particles  will  leave  the  active 
region  of  the  cloud  earlier  due  to  centrifugal  or  gravitational 
foxes.  But  because  the  model  divides  the  fission  product  Inventory 


A-2 


Into  a  surface  distributed  component  and  a  volume  distributed 
component.  It  can  be  used  to  develop  a  number  of  composite  activity 
size  distributions  for  the  AFIT  model.  More  specifically.  If  S(d)  Is 
the  distribution  of  activity  on  the  surfaces  of  the  particles  and  Y(d) 
is  the  distribution  of  activity  inside  of  the  particles,  and  If  the 
particle  size  distribution  N(d)  is  taken  to  be  lognormal,  then  the 
total  activity  size  distribution  takes  on  a  vei7  powerful  form  because 
of  the  properties  of  the  lognormal  distribution.  Surface  area  Is 
distributed  as  the  second  moment  of  the  particle  size  distribution  and 
volume  is  distributed  as  the  third  moment  of  the  particle  size 
distribution.  Thus  the  total  activity  size  distribution  becomes 

A{d)  .  Aj  S{d)  ♦  A^  V(d)  (A-1) 

or 

A(d)  «  Aj  LNCag.B)  ♦  A^  LNCaj.s)  (A-2) 

where 

Aj  «  the  surface  distributed  activity 

Ay  a  the  volume  distributed  activity 

LN(aj,B)  a  the  lognormal  distribution  with  median  aj  and 


standard  deviation  b 


Is  the  median  for  the  n-th  moment  of  the  particle  size  distribution. 
Since  the  two  parts  of  the  distribution  have  their  activity  computed 
separately,  this  allows  for  two  K-factor  constants  (these  are  used  to 
Convert  from  activity  to  dose  rate)  to  be  computed  at  the  time  of 
Interest.  The  effect  of  this  Is  to  allow  for  differences  In  the 
average  energy  of  emission  for  the  two  groups.  In  addition,  the  decay 
Is  usually  computed  with  the  Way-WIgner  formula  which  Is  approximately 
correct  for  fission  products  taken  as  a  whole,  but  when  the  debris  Is 
fractionated,  this  method  Is  no  longer  applicable.  Prelling  has  shown 
that  at  late  time  there  Is  considerable  difference  between  the  rates 
of  decay  for  the  two  groups  (Frelling,  1964:  7-10).  This  formulation 
allows  for  separate  approximate  rates  of  decay  (perhaps  curve  fits 
from  a  parametric  study)  to  be  used  for  each  group. 
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R|  j  Plots  for  Small  Boy 

This  appendix  Includes  fractionation  plots  for  all  chains  for 
which  data  Is  available  In  the  literature  (Crocker,  196S:  72-81). 
nuclides  used  for  radiochemical  analysis  are  as  follows: 
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Least  squares  fits  to  the  test  data  are  shown  on  the  plots  as  the. 
solid  straight  lines.  Mass  chains  89  and  95  are  reference  chains, 
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Ffgure  23.  Snail  Boy  Fractionation  Plots:  Mass  Chain  137 
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Figure  27.  Snail  Boy  Fractionation  Plots:  Mass  Chain  136 


Figure  28.  Small  Boy  Fractionation  Plots:  Mass  Chain  140 
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j  Plots  for  Johnny  Boy 

This  appendix  Includes  fractionation  plots  for  all  chains  for 
which  data  Is  available  In  the  literature  (Crocker,  1965:  72-81).  The 
nuclides  used  for  radiochemical  analysis  are  as  follows: 
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Least  squares  fits  to  the  test  data  are  shown  on  the  plots  as  solid 
straight  lines.  Mass  chains  89  and  95  are  reference  chains. 
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Specific  Activity  Plots  for  Small  Boy 

This  appendix  Includes  plots  of  the  specific  activity  for  all 
chains  of  general  Interest.  Mass  chains  which  are  volume  distributed 
will  have  constant  specific  activity  curves.  Mass  chains  which  are 
surface  distributed  will  show  a  marked  decrease  In  specific  activity 
as  particle  size  Increases.  Chains  which  have  mixed  behavior  will 
exhibit  a  combination  of  these. 
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Figure  47.  Snail  Boy  Specific  Activity  for  Mass  Cham  132 
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Figure  55.  Small  Boy  Specific  Activity  for  Mass  Chain  144 
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